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DISPOSITIF ELECTRONIQUE DE TRAITEMENT DE DONNEES-IMAGE, POUR LA SIMULATION DU 
COMPORTEMENT DEFORMABLE D'UN OBJET. 

Un dispositif 6lectronique comprend una interface-uti- 
lisateur a retour de force (4) conformement aux reactions 
d'un outil, un module de "collision" (18) propre a estimer un 
point d' intersection entre une droite materialisant un depla- 
cement tire d'une action de I'outil et le maillage surfacique 
d'un objet designe, et un module de forces internes (16) qui 
etablit un champ de forces Internes entre des noeuds d'une 
partie au moins du maillage volumique de I'objet a partir du 
deplacement applique aux noeuds appartenant a la mallle 
surfacique contenant le point d'intersection, de conditions 
aux llmites, et de tenseurs de noeud et de tenseurs de 
liaison issus respectivement pour chaque noeud et chaque 
liaison de cette partie au moins, de matrices de rigidity, et 
un module de reaction (20) pour determiner la force de reac- 
tion de I'objet correspondant h sa deformation estim^e a 
partir des forces internes, de sorte que la force generee par 
rinterface-utilisateur (4) soit equilibree par cette force de 
reaction. 




STATION OUPHOJE 



uNmoEOOKnm 



BNSDOCID: <FR 2771202A1_L> 



2771202 



n-ic pr>gHiif ^lecuroniaue de tra^^emen i- rie donn6e s- image . pour 14 
cimnlanion du comportement: d^formable d'un obiet. 

5 

L' invention concerne le domaine du traitement de donn^es- image 
n\am6riques d'un jeu repr6sentatif d'une image tridimensionnelle 
(3D), pour la simulation du comportement dSformable d'un objet. 

10 L' invention s' applique plus particulierement , mais de fa9on non 
exclusive, au traitement d'un jeu de donn^es-image d'une image 
dite m^dicale. 

Dans de nombreux domaines, il est trfes int6ressant de pouvoir 
15 simuler des interventions d'un op6rateur, & I'aide d'un (ou 
plusieurs) outil(s) connu(s), sur un (ou plusieurs) objet(s) 
d6f ormable ( s ) . On entend ici par intervention, aussi bien une 
manipulation, en vue par example d'un d6placement, qu'une 
transformation locale, comme par exemple, dans le cas d'une 
20 intervention chirurgicale, de type incision ou extraction d'une 
partie d'un organe. 

La simulation consiste h. afficher 1' image d'un objet, et 
6ventuellement de la region dans laquelle il se trouve habituel- 

25 lement, et la representation d'un outil au moins dont le 
d6placement "virtuel", relativement 1' objet, est d6fini par 
une interface-utilisateur dont un ensemble-harnachement est 
manoeuvre par un op6rateur, en vue de la simulation du maniement 
dudit outil. Afin de pouvoir simuler la reaction de 1' objet sur 

30 1' outil, 1 ' interface-utilisateur est capable de gen6rer un 
retour de force, conformSment aux reactions de 1' outil. On 
entend par force de reaction d'un objet le retour de force. 

Dans les dispositifs connus, un module de reaction permet de 
35 determiner cette force de reaction de 1' objet h. partir d'une 
deformation estimee de cet objet. Cette deformation s'obtient h 
I'aide d'un module de forces internes et d'un module de 
rafraichissement d' image. Le module de forces internes est 
capable, sur designation d'un objet 3D apparaissant dans un jeu 
40 de donnees-image, d'6tablir ian champ de forces internes. 
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repr6sentatif de la deformation de I'objet, entre des noeuds 
d'un maillage volximique fonction d'un maillage surfacique de cet 
Ob jet, ^ partir d'une loi de deformation et d'une action d^finie 
par 1' interf ace-utilisateur et representative d'une manoeuvre de 
5 I'outil. 

Le module de raf raichissement permet alors de calculer de 
nouvelles donn^es-image de I'objet^ en presence des deformations 
estimees, compl6t6es de la representation de I'outil. Ces 

10 nouvelles donnees-image qui forment la nouvelle image de 
I'objet, et eventuellement celle de la region qui I'entoure, 
sont alors affichees sur un dispositif de visualisation, de 
sorte que I'operateur puisse voir en temps reel le resultat de 
la manipulation de 1 ' ensemble-harnachement qui simule 1' action 

15 sur I'outil. 

Un tel dispositif doit permettre la formation d'un operateur ou 
bien la mise au point de nouvelles techniques d' intervention sur 
1' Ob jet. Dans un domaine tel que la chirurgie, et plus particu- 

20 lidrement encore dans le domaine de la chirurgie laparoscopique , 
ce type de dispositif peut permettre de sauver des vies 
humaines. Pour ce faire, il est imperatif que la simulation 
permette de reproduire le plus f idfelement possible le geste de 
1' operateur (ou en d'autres termes son action sur un outil, ici 

25 virtuel). Cela necessite un traitement en temps reel des 
donnees-image couple a une restitution des forces induites par 
I'objet en reponse ^ la deformation gen6ree par "I'outil". 

Or, du fait des techniques de calcul utilisees par les disposi- 
30 tifs connus, 1' estimation des forces internes n6cessite des 
temps de calcul importants incompatibles avec une simulation 
dynamique continue^ En d'autres termes, les dispositif s actuels 
ne permettent pas de visualiser, de faqon continue pour un oeil 
humain, la totalite de 1' action d'un outil sur un objet 
35 deformable. 
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Par ailleurs, aucun dispositif actuel ne permet de simuler en 
temps r6el une action de type incision, ou d6chirement , ou bien 
encore retrait de matiere d'un objet d6formable. 

5 La pr6sente invention a done pour but de r6soudre tout ou partie 
des inconv6nients pr6cit6s dans le domaine du traitement de 
donn^es -image numSriques d' objet 3D. 

Elle propose k cet effet un dispositif 61ectronique de traite- 
10 ment de donn6es-image du type d6crit dans 1 ' introduction, dans 
lequel, d'une part, on pr6voit un module de "collision" capable 
d'estimer un point d' intersection entre une droite mat^rialisant 
un d^placement tir6 de 1' action definie et le maillage surfaci- 
que, et d' autre part, le module de forces internes est agence 
15 pour estimer la force interne exerc6e sur chaque noeud d'une 
premiere partie au moins du maillage volumique de 1' objet ci 
partir du d^placement tir6 de 1' action, appliqu6 aux noeuds 
appartenant St la maille surfacique contenant le point d' inter- 
section^ de conditions aux limites^ et de tenseurs de noeud et 
20 de tenseurs de liaison issus respect ivement pour chaque noeud et 
chaque liaison d'au moins la premiere partie au moins, de 
matrices de rigidity propres a chaque maille volumique d'au 
moins la premiere partie et fonction de la loi de deformation. 

25 Bien entendu, la premiere partie de la maille volumique h 
laquelle est appliquee la technique ci-dessus, que I'on 
appellera par la suite "masses/tenseurs" , pourra etre 6gale a la 
maille volumique complete. Dans le cas contraire (lorsqu'il 
s'agit ef f ectivement d'une partie de cette maille volumique), 

30 les forces internes appliqu6es aux noeuds de la partie compl6- 
mentaire de cette premiere partie (appel^e par exemple seconde 
partie) sont d6termin6es k partir d'une autre technique, comme 
par exemple celle des 616ments finis reposant sur des pr6- 
calculs memorises de fagon Si permettre des calculs en tenths 

35 r6el. Une telle technique est enseign6e notamment dans 1' article 
de S. Cotin, H. Delingette, M. Bro-Nielsen et N. Ayache, 
"Geometric and physical representations for a simulator of 
hepatic surgery", public dans les actes de la Conference 
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Medecine meets with virtual reality de janvier 1996. Dans ce qui 
suit, on appellera module hybride la double technique de calcul 
des forces internes et de la deformation de I'objet. 

5 Cette technique dite des masses/tenseurs utilisSe pour le calcul 
des forces internes et de la deformation de la maille volumique 
de I'objet autorise la simulation en continu d'une action 
exerc^e sur un outil virtuel au moins . II est clair que le temps 
de calcul sera d'autant plus r6duit que la dimension de la 
10 premiere partie de la maille volumique sera petite. 

Selon une autre caract^ristique de 1' invention^ le dispositif 
peut comprendre un module de maillage lui permettant de designer 
le ou les objet(s) 3D sur lesquels la simulation doit etre 

15 effectu6e par une determination d'enveloppe externe, puis de 
decomposer cette ou ces enveloppe(s) en mailles surfaciques, de 
preference de forme triangulaire, et enfin de decomposer le 
volume interne de chaque enveloppe en mailles volumique & partir 
du maillage ^urfacique correspondant de ^agon fournir le 

20 maillage volumique de I'objet associe. II est clair que dans le 
cas prefere d'un maillage surfacique triangulaire les mailles 
volimiques seront de forme tetraedrique. Ces formes sont 
actuellement preferees car elles permettent une modeiisation 
precise d'un objet de forme complexe. Mais bien entendu, 

25 d'autres type de maillage peuvent etre utilises. 

De telles enveloppes externes et mailles volumiques peuvent etre 
obtenues par des methodes respectivement de segmentation (par 
exemple par extraction d'iso-surfaces) et de type Delaunay- 
30 Voronoi, Toutes ces methodes sont bien connues de I'homme de 
I'art. 

Dans un mode de realisation du dispositif, son module de forces 
internes est capable de calculer lui-meme les matrices de 
35 rigidite de chaque maille volumique, ainsi que les tenseurs de 
noeud et tenseurs de liaison. 
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Ces calculs sont, coinine cela a 6t6 dit pr6c6cieininent , effectues 
k partir d'une loi de deformation qui est de pr^f^rence, de type 
61astique lineaire volumique. En d'autres termes^ la force 
exerc6e sur un noeud depend des d6placements respectivement de 
5 ce noeud et des noeuds auxquels il est connect^, relativement h 
leurs positions d'6quilibre respectives- Bien entendu, on 
pourrait utiliser d'autres lois de deformation plus complexes et 
notaniment non-lineaires . 

10 Le module de forces internes pourra 6galement §tre agenc6 pour 
determiner les forces internes exerc6es sur certains au moins 
des noeuds de la premi&re partie du maillage volumique h partir 
de la loi de deformation et de forces auxiliaires surfaciques 
dependant de paramStres choisis de I'objet, stock6s, comme par 

15 exemple la texture de I'objet, la presence de structures et 
sous-structures sous- jacentes, etc* 

On entend ici par forces auxiliaires surfaciques^ par exemple 
des tensions de surface qui perroettront^ dans certaines 
20 situations telles qu'une incision, d'amplifier un effet visuel 
au niveau de I'affichage. 

De meme, le module de forces internes pourra etre agenc6 pour 
estimer les d6placements des noeuds du maillage volumique (au 

25 moins sa premiere partie) a partir du d6placement tire de 
1' action d6finie et de forces externes, en particulier de type 
force gravitationnelle et/ou forces d' interaction entre objets 
d'une meme region. Cela permet de tenir compte, d'une part, de 
I'af f aissement partiel d'un objet sous son propre poids, et 

30 d' autre part, de la presence d' objets voisins et des liens qui 
existent avec ces objets voisins - 

Preferentiellement, les d6placements estim6s des noeuds, autres 
que ceux de ladite maille surfacique comprenant ledit point 
35 d' intersection, sont calcul6s a partir des forces internes par 
integrations successives & I'aide d'une methode choisie parmi au 
moins la methode d'Euler et la methode de Runge-Kutta, et plus 
preferentiellement encore par la methode de Runge-Kutta dite 
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"d'ordre 4". Bien entendu, d'autres methodes d ' integration 
peuvent etre envisagSes. 

Selon encore une autre caracteristique de 1' invention , le module 
5 de forces internes peut etre capable de simuler des deformations 
non seulement de type g6om6trique, mais ^galement de type 
incision et/ou retrait de matifere et/ou d6chirement. 

Pour permettre la simulation de d6coupe (ou incision) et/ou de 
10 d6chirement (ou fracture), le module de forces internes est 
propre, apr^s determination des d^placements estim^s des noeuds, 
ci supprimer au moins une liaison entre des noeuds voisins en 
fonction d'un premier critere, puis a mettre & jour les tenseurs 
de noeud et les tenseurs de liaisons en fonction de la ou des 
15 liaisons supprimees^ et enfin h recalculer les forces internes 
des noeuds d'au moins la premifere partie du maillage volumique. 

De preference, le premier critere porte sur au moins un 
parametre choisi parmi au moins une information transmi.se par 
20 I'interf ace-utilisateur et relative au type d'outil manoeuvre, 
une variation de volume de la maille volumique comprenant la 
liaison h supprimer, et une variation de longueur d'une liaison 
de la maille voluiaique comprenant la liaison h supprimer • 

25 On entend ici par information, par exemple une donnee precisant 
que I'outil est manoeuvre en vue d'une incision ou d'une 
destruction de mati&re. 

De meme, pour permettre la simulation de retrait de matiere, le 
30 module de forces internes est propre, apr^s determination des 
deplacements estimes des noeuds, h supprimer un noeud en cas de 
detection de suppression de toutes les liaisons qui relient 
ledit noeud aux noeuds voisins ou en fonction du premier 
critere, puis h mettre ci jour les tenseurs de noeud et les 
35 tenseurs de liaisons en fonction du noeud et des liaisons 
supprimes, et enfin ^ recalculer les forces internes des noeuds 
d'au moins la premiere partie du maillage volumique. 
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Les outils (ici virtuels) capables de d^coupe (ou d' incision) 
et/ou de retrait de mati^re sont par example des scalpels, 
pinces coupantes, ou bien bistouri m6canique ou 61ectrique, ou 
bien encore lasers. 

5 

De plus, le module de forces internes est de pr6f6rence agenc6 
pour, en cas de suppression de liaison et/ou de noeud et avant 
mise & jour des tenseurs de liaison et de noeud, a j outer de 
nouveaux noeuds ind^pendants entre eux et de nouvelles liaisons 
10 de maniere H remailler localement le maillage volumique suite a 
la suppression. 

Lorsque le dispositif est agenc6 pour travailler selon le modfele 
hybride pr6cit6, son module de forces internes est capable de 

15 determiner les forces internes exerc6es sur les noeuds d'au 
moins une seconde partie du maillage volumique h partir de 
conditions aux limites d^finies par des noeuds dits "de 
connexion" places h 1' interface entre les premifere et seconde 
parties, et d'une table de tenseurs de deformation, dont chaque 

20 tenseur est repr6sentatif de 1' influence d'un d6placement 
eiementaire de chaque noeud d'au moins la seconde partie sur 
chaque autre noeud d'au moins cette seconde partie. 

Les conditions aux limites servant au calcul des forces internes 
25 de la seconde partie sont pr6f6rentiellement definies par les 
forces internes calcul6es pour les noeuds de connexion lors du 
calcul des forces internes des noeuds de la premiere partie. 

Dans le module hybride, le module de forces internes procSde de 
30 preference par iterations successives jusqu'i obtention d'une 
position dite "d'equilibre- des forces internes des noeuds de 
connexion. Pour ce faire, ce module est agence pour deduire des 
valeurs des forces internes exerc6es sur les noeuds de la 
seconde partie du maillage volumique, des valeurs de deplacement 
35 des noeuds de connexion de manifere h fournir des conditions aux 
limites permettant ^ leur tour de calculer les forces internes 
des noeuds de la premifere partie. 
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La subdivision du maillage volumique en parties (au moins deux) 
est d^termin^e ^ partir d'un critere pr6d6termin6 portant au 
moins sur un parametre des donnees- image de I'objet choisi parmi 
des param^tres physiques et des paramStres anatomiques, en 
5 particulier sur une intensite. Cette subdivision peut Stre 
effectu6e pendant la formation du maillage volumique de I'objet 
3D- Par consequent, lorsque le dispositif ne comprend pas de 
module de maillage, la subdivision peut etre effectu6e par un 
traitement externe. Le dispositif peut comprendre un module de 
10 partitionnement , par exemple faisant partie du module de 
maillage si celui existe, destin6 2i fournir la subdivision du 
maillage volumique. 

Pr6f6rentiellement, le module de determination de "collision" 
15 est agence pour determiner une collision entre au moins deux 
outils geres par 1 ' inter face-utilisateur- Ce qui permet de g6rer 
des conflits lorsque un ou plusieurs operateur(s) manoeuvre ( nt ) 
au moins deux outils en meme . temps. 

20 Dans ce but, chaque outil est represents par au moins un point 
materialisant son extr6mite d' interaction avec I'objet, et une 
multiplicite de points relies entre eux, ainsi qu'a I'extremite, 
par des segments, materialisant sa "tige". 

25 De preference, le module de determination de collision estime 
les coordonnees du point d' intersection entre "1' outil" et une 
maille surfacique comme suit. 

Tout d'abord, il cree un espace tridimensionnel englobant 
30 I'enveloppe externe de I'objet, puis il decompose cet espace en 
cellules volumiques dont le nombre est choisi de sorte que 
chaque cellule comporte un nombre de noeud du maillage volumique 
de I'objet sensiblement egal au nombre de noeuds contenus dans 
les autres cellules volumiques, chaque cellule intersectant la 
35 surface externe comportant au moins un noeud, et enfin il stocke 
dans des multiplets les coordonnees de chaque noeud. en reference 
k la cellule volumique qui I'englobe. La presence d'un point de 
1' outil dans 1' espace s'effectue alors, de preference, par une 
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comparaison entre les multiplets et les coordonn^es du point. 
L'un de ces multiplets permet de designer la cellule volumique 
de I'espace dans laquelle se trouve le point. 

Puis^ de preference ^ le module de detection de collision 
determine la distance (euclidienne) qui s6pare le point de 
I'outil du ou des noeuds englob^s dans la cellule volumique 
designee de fagon h determiner la plus petite de ces distances, 
appeiee distance minimale. II determine ensuite la distance (par 
exemple euclidienne) qui s6pare le point du ou des noeuds 
eng lobes dans un nombre predetermine de cellules volumiques 
voisines de la cellule volumique dans laquelle il se trouve de 
fagon h comparer ses distances k la distance lainimale. Puis, il 
determine 1' ensemble des mailles surfaciques adjacentes au noeud 
associe ^ la distance minimale de fagon a determiner si un 
segment defini par la position du point de I'outil et par sa 
position precedente intersecte I'une au moins de ces mailles 
surfaciques adjacentes. Il en d6duit finalement les coordonnees 
(barycentri^jues) du point d ' intersection (ou de collision) entre 
1' Ob jet et I'outil, en vue de leur transmission au module de 
forces internes - 

De fagon preferentielle, le module de detection de collision est 
capable de modifier le contenu des multiplets entre deux 
25 determinations de presence de points de I'outil h I'interieur 
des cellules volumiques, en cas de detection d'une deformation 
du maillage par le module de forces internes. Cela permet 
d'ameiiorer la precision de la detection de collision et done du 
calcul de la deformation. 

30 

Afin de permettre une simulation d' action sur un ou plusieurs 
outils la plus proche de la r^alite, 1 ' interf ace-utilisateur 
comprend un ensemble-harnachement manoeuvrable par au moins une 
main d'operateur pour simuler la manoeuvre de chaque outil. 
35 L'op6rateur manoeuvre ainsi 1' ensemble-harnachement^ qui peut 
etre un "joystick- muni 6ventuellement d ' actionneur ( s ) ou bien 
un "gant" articul6 fixe sur sa main, puis 1 ' interf ace-utilisa- 



10 



15 



20 
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teur d^finit le deplacement de L'outil virtuel associ^ h partir 
de cette manoeuvre. 

Salon encore une autre caract^ristique de 1' invention, le 
5 dispositif comprend des moyens d'affichage permettant de 
visualiser en temps r6el (et en continu) 1' image ( f orm6e des 
donn^es- image) de I'objet et d'une representation de I'outil. 

L' invention s' applique tout particuliferement aux jeux de 
10 donn^es*- image reprSsentant une image num^rigue tridimensionnelle 
d'une region comportant au moins un objet 3D, dont I'objet 
d^sign^, et plus particuliferement encore aux jeux de donnSes- 
image repr^sentant une image numSrique tridimensionnelle d'une 
region d'un etre vivant (animal ou humain) comprenant des 
15 structures anatomiques deformables telles que le foie, le rein, 
la v6sicule biliaire, ou bien encore le coeur. 

L' invention propose egalement un proc6d6 de traitement de 
donn^es— image num6ric[ues pour la mise en oeuvre du dispositif 
20 dScrit ci-avant, comprenant les Stapes connues suivantes : 

- pr6voir une interf ace-utilisateur capable de g6n6rer un retour 
de force, conform^ment aux reactions d'un outil, 

- 6tablir, ^ partir d'une loi de deformation et d'une action 
definie par 1' interf ace-utilisateur et representative d'une 

25 manoeuvre de 1' outil, un champ de forces internes entre des 
noeuds d'un maillage volumique fonction d'un maillage surfacique 
d'un objet 3D apparaissant dans un jeu de donnees-image, 

- determiner la force de reaction de I'objet qui correspond a sa 
deformation estimee h partir des forces internes, de sorte que 

30 la force gen6ree par 1 ' interf ace-utilisateur soit sensiblement 
equilibree par cette force de reaction, 

- calculer de nouvelles donnees-image de I'objet, en presence 
des deformations estimees, compietees de la representation dudit 
outil, 

35 et caracterise par le fait que I'on pr6voit une etape dans 

laquelle on estime un point d ' intersection entre une droite 
materialisant un deplacement tire de 1' action definie et le 
maillage surfacique, et que I'on estime la force interne exerc6e 
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sur les noeuds d'une premiere partie au moins du maillage 
volumique de I'objet, ci partir du deplacement applique aux 
noeuds appartenant h la maille surfacique contenant le point 
d' intersection, de conditions aux liitiites , et de tenseurs de 
5 noeud et de tenseurs de liaison issus respectivement pour chaque 
noeud et chaque liaison de cette partie au moins, de matrices de 
rigidity propres h chaque maille volumique d'au moins la 
premifere partie et fonction de la loi de deformation. 

10 D'autres caract6ristiques et avantages de 1' invention apparal- 
tront h I'examen de la description detaill6e qui va suivre, 
ainsi que des dessins annexes , sur lesquels : 

- la figure 1 est une photographie d'une partie d'un dispositif 
15 de simulation selon 1' invention; 

- la figure 2 est un schema fonctionnel illustrant 1' architec- 
ture d'un dispositif de simulation selon 1' invention dans un 
mode de r6alisai:ion h deux modules de -traitemeni: ; 

20 

- la figure 3 illus-tre sch^matiquement le mode de deplacement 
d'un outil pilots par un type d' inter face-utilisateur du 
dispositif selon 1' invention; 

25 - la figure 4 illustre I'enveloppe externe (ou surface externe) 
d'un foie humain obtenue ci partir d'une technique de segmenta- 
tion; 

- la figure 5 illustre I'enveloppe externe du foie de la figure 
30 4 muni d'un maillage surfacique de type triangulaire; 

- la figure 6 est un schema illustrant un maillage volumique 
t6tra6drique obtenu partir d'un maillage surfacique triangu- 
laire du type de celui de la figure 5; 

35 

- la figure 7 est un schema illustrant les d6placements h un 
instant t de deux noeuds par rapport k leurs positions d'6quili- 
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bre respectives et des positions d'6quilibre de leurs noeuds 
voisins; 

- la figure 8 est un schema illustrant une maille volumique 
5 t6traedrique et la designation des tenseurs de noeud et des 

tenseurs de liaison respectivement associ^s ^ ses noeuds et 
liaisons; 

- la figure 9 est un exemple de matrice de rigidity d'un objet 
10 3D repr6sent6 par 15 noeuds, les croix (X) signalant des valeurs 

non nulles et les points (.) des valeurs nulles; 

- la figure 10 est un m6canisme illustrant les Stapes de 
determination d'un maillage final d'un objet 3D ^ l'6quilibre; 

15 

- la figure 11 est un m^canisme illustrant les Stapes de 
determination d'une deformation et de 1' image de 1' objet 
def orme ; 

20 - la figure 12 est un schema illustrant la modification de la 
topologie d'un maillage volumique tetraedrique suite h la 
rupture d'une liaison; 

- les figures 13A a 13D illustrent une simulation sur quatre 
25 images successives de 1' action de l'extr6mite d'un outil virtuel 

sur un lobe de foie; 

- la figure 14 est un schema illustrant une decomposition d'un 
maillage surf acique triangulaire en deux domaines (ou parties), 

30 dans un traitement de type hybride; 

- la figure 15 est un mecanisme illustrant une boucle de 
determination de positions d'equilibre dans un maillage hybride; 

35 - la figure 16 est un schema illustrant un espace englobant le 
foie de la figure 4 et decompose en cellules volumigues 
paralieiepipediques ; 
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- la figure 17 est un schema illustrant une partie d'un maillage 
volumique t6tra6drique contenu dans une cellule volumique 
parall616pipedique, ainsi qu'un point de collision entre un 
outil et I'une des mailles surfaciques triangulaires de la 
partie dudit maillage volumique; et 

- la figure 18 est un schema illustrant la materialisation d'un 
outil sous forme de points et segments, en reference h. une 
partie d'un maillage surfacique triangulaire - 



10 



Les dessins sont, pour I'essentiel, de caractfere certain. En 
consequence, ils font partie int^grante de la pr6sente descrip- 
tion, lis pourront done non seulement servir h mieux faire 
coitiprendre 1' invention, mais aussi contribuer h la definition de 
15 celle-ci. 

Dans ce qui suit, on fera reference Si un traitement d'un jeu de 
donnees-image num^riques formant une image tridimensionnelle 
(3D) medicale, et plus particuli^rement , mais seulement A titre 
20 d'exemple, h des images de regions du foie du type de celle 
illustree partiellement sur les figures 1 et 4, lesquelles ont 
6t# obtenues sur un sujet humain. 

II est cependant clair que 1' invention n'est pas limitee au 
25 traitement d'images m^dicales, et encore moins h celui d'images 
du foie. Elle s' applique de fagon g^nerale au traitement 
d'images numfiriques d'objets tridimensionnels (3D), notamment 
deformables, en vue d'une simulation de leur deformation par au 
moins un outil virtuel. Le mot deformation est ici entendu dans 
30 son acceptation la plus generale, h savoir aussi bien de type 
geometrique (surfacique ou volumique) que de type incision, ou 
dechirement, ou bien encore retrait de matiSre. 

En matifere m6dicale, un jeu de donn6es-image formant une image 
35 tridimensionnelle (3D) est obtenu a I'aide d'appareils tels que 
les scanners h rayons X ou les appareils de resonance magnetique 
^ nucieaire (IRM), ou plus g6neralement tout type d'appareil 
capable d'acquerir des images i variation d'intensite. Chaque 
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partie 616mentaire d'une region representee par une image 3D est 
d6finie par 3 coordonn^es spatiales et au moins une grandeur 
physique mesur6e, comme par exemple une intensity. 

5 Dans le cas d'une IRM^ 1' image 3D d'une region observee est 
constituee d'une multiplicity de coupes 2D eitpil6es, dans 
lesquelles les variations d' intensity repr6sentent la density en 
proton des tissus. 

10 Sur la figure 1 se trouve partiellement illustr6 un dispositif 
selon 1' invention destine ^ simuler le comportement dynamique 
d'un objet tridimensionnel (3D). Ici, comme indiqu6 prec^dem- 
ment, 1' objet 3D est un foie humain, 

15 Le dispositif comprend dans 1' exemple illustre, une unite de 
controle et de calcul 1 coupiee une station graphique 2 
comprenant des moyens d'affichage, comme par exemple un moniteur 
3, ainsi qu*k une interf ace-utilisateur 4 Sl retour de force. 

20 Dans cet exemple^ 1' unite de contraie 1 et la station graphique 
2 sont deux calculateurs s^pares, mais il est clair qu'ils 
pourraient former une unique unite de calcul. 

Dans le mode de realisation illustre, 1 'unite de controle 1 est 
25 de preference un micro-ordinateur de type PC, tandis que la 
station graphique 2 est de preference une station de travail de 
type DEC ALPHA 233 MHz 6quipee d'une carte graphique 3D 
acceieratr ice . Il est ici preferable d'utiliser 1' unite de 
calcul la plus puissante pour gerer I'affichage des images, du 
30 fait de la n6cessite de calculer la deformation et d'afficher en 
temps reel ^ la fois 1' image du foie et de la representation 
d'un outil, ou plusieurs, ce qui necessite un raf ralchissement 
des images & une frequence de type video super ieure & 20 Hz. 

35 Par ailleurs, la liaison entre I'unite de contrdle 1 et 
1' interf ace-utilisateur Si retour de force 4 sera de preference 
un bus ISA, lequel autorise des vitesses de transfert de donnees 
de 10 kHz bien superieures ^ la frequence requise pour le retour 
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de force, laquelle est g6n6ralement comprise entre 300 et 500 
Hz. 

L' inter face-utilisateur a retour d' effort 4 comprend une 
coinmande manuelle 5 (ou ensemble-harnachement ) de type "joys- 
tick" capable de reproduire le comportement d'un outil de 
chirurgie minimalement invasif . II pourra s'agir notamment d'un 
outil de chirurgie laparoscopique ou coelioscopique- Dans ce qui 
suit/ on consid6rera que le joystick simule le mouvement d'une 
pince coupante microscopique destin6e h etre introduite dans le 
corps d'un patient via un point fixe appele trocart en vue d'une 
intervention sur le foie. Bien entendu, 1' invention n'est pas 
limitee & ce seul type de chirurgie, ni d' outil. 

L' inter face-utilisateur §i retour d' effort 4 pourra par exemple 
etre le "laparoscopic impulse engine" (LIE) de la Soci^te 
Immersion Corporation, ou bien le dispositif PHANTOM d6velopp6 
par le Massachussets Institute of Technology. Le "laparoscopic 
impulse engine" (LIE) possede cinq degree de liberty, trois pour 
repr^senter les d6placements dans I'espace de la tige de 
1' outil, un quatrifeme pour la rotation de cette tige autour d'un 
axe passant par le trocart, et un cinquifeme associ^ a I'ouver- 
ture ou la fermeture de la pince plac6e k I'extr6mit6 de la 
tige. Le retour d' effort peut etre appliqu6 Si trois de ces 
degr6s de liberty h I'aide de moteurs. La transmission des 
forces est r6alis6e via un systfeme de cables et de cabestans 
d61ivrant ainsi un couple elev6 et progressif- 

Sur la figure 3 se trouve illustr^ sch^matiquement le mode de 
30 d6placement de la tige 6 de la pince coupante pilot6e par le 
LIE- L'extr6mit6 7 de la tige 6 est ici mat6rialis6e par une 
sphere qui repr^sente schema tiquement I'extr6mit6 active (par 
exemple coupante) de 1' outil virtuel. Le trocart est materialist 
par un point fixe 8 par lequel passe la tige 6 de 1' outil, 
35 laquelle est connectte h deux moteurs 9 et 10 qui lui permettent 
des d6placements sur une calotte sph6rique dont le rayon varie 
selon la longueur de la tige 6, laquelle est g6r6e par un 
troisifeme moteur non repr6sent6. La position de I'extr6mit6 7 de 
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I'outil sur la calotte sph^rique est d6tenriin6e par les valeurs 
des angles a Bt, fi. 

Lorsqu'un op6rateur manoeuvre, k I'aide de sa main 11 (voir 
5 figure 1), la partie pr6hensible du joystick 5, il peut d^placer 
I'axe 6 de mani^re h faire varier la longueur qui s6pare le 
point fixe 8 (ou trocart) de I'extr6mit6 coupante 7 permettant 
ainsi de simuler le rapprochement de I'outil vers I'objet 3D, 
mais il peut 6galement dSplacer cette extr6mit6 7 sur une 
10 calotte sph6rique dont le rayon depend de la distance trocart- 
extr6mit6 • 

Bien entendu, 1 ' interf ace-utilisateur 4 peut comprendre deux LIE 
de fagon Si permettre la simulation du maniement simultan6 de 
15 deux outils- 

Un but du dispositif selon 1' invention est, coxnme cela a d^j^ 
6t6 dit, de traiter en temps r6el des donnSes-image n\am6riques 
formant des images successive^ tridimensionnelles d'un objet SD 
20 (ici un foie humain), obtenues, par exemple, par scannograpie 
par rayons X, ou par IRM, selon I'objet ou organe concern^. 

II faut tout d'abord que le dispositif possdde le jeu de 
donn6es-image foirmant 1 ' image de I'objet 3D, I'enveloppe externe 

25 de cet objet, le maillage surfacique de cette enveloppe, et le 
maillage volumique primaire de I'objet form6 ^ partir du 
maillage surfacique. Toutes ces donn^es peuvent etre obtenues en 
dehors du dispositif selon 1' invention. Elles pourront notamment 
provenir d'une banque de donn^es d' images coupl^e ^ un module de 

30 maillage. Cependant, le dispositif peut comprendre un module de 
maillage capable de determiner 1 'enveloppe externe et les 
maillages surfacique et volumique primaire sur la base d'un jeu 
de donn^es- image fourni par un moyen externe. 

35 L'enveloppe externe 12 (voir figure 4) pourra fitre obtenue par 
toute technicjue connue de I'homme de I'art, et par exemple i 
I'aide d'une technique de segmentation de type extraction d'iso- 
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surface telle que les techniques dites des "SD snakes" ou bien 
des "marching cubes", bien connues de I'hoirane de I'art. 

Dans le cas d'un module de maillage integre au dispositif, 
5 celui-ci sera de pr6f6rence log6 dans la station graphique 2. 

La forme g6n6rale des mailles surfaciques 13 d6composant I'enve- 
loppe externe est de pr6f6rence triangulaire (voir figure 5). Un 
critfere d' optimisation de la surface des mailles surfaciques 

10 pourra §tre mis en oeuvre afin d'obtenir un maillage r^gulier, 
Bien entendu, d'autres types de maillages surfaciques peuvent 
etre utilises, corome notamment les maillages "simplexes" 
d6velopp6s par H, Delingette et divulgues dans le rapport 
technique de I'INRIA N*" 2214 sous le titre "Simplex meshes : a 

15 general representation for 3D shape reconstruction", 1994. 

A partir des mailles surfaciques 13, on determine un maillage 
volumique de I'objet 3D, ici le foie. II est clair que la forme 
g^n^rale des mailles depend de la forme des maille^s surfaoi<iues 
20 13. Dans cet exemple ou les mailles surfaciques 13 sont 
triangulaires , les mailles volumiques (voir figure 6) sont 
avantageusement des t6tra6dres 14 aux quatre sommets desquels se 
trouvent places des noeuds 15 sur lesquels est r6partie la masse 
de I'objet, laquelle est connue par ailleurs. 

25 

Cette decomposition du volume englob^ par I'enveloppe externe 12 
en mailles, ici des t^trafedres, peut etre obtenue h I'aide d'un 
algorithme de type Delaunay-Voronoi utilisant, par exemple, le 
module Simail commercialism par la Soci6t6 frangaise SIMULOG 
30 S.A. Bien entendu, d'autres algorithmes pourront etre utilises 
en vue de I'obtention du maillage volumique. 

Comme cela a dSj^ 6t6 indiqu6 pr6c6demment , le dispositif salon 
1' invention ne coitprend pas n6cessairement un module de maillage 
35 destine & determiner I'enveloppe externe de I'objet 3D et/ou la 
decomposition en mailles surfaciques et/ou la decomposition en 
mailles volumiques. Tout ou partie de ce traitement du jeu de 
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donn^es- image de I'objet 3D peut etre affectum de fagon 
totalement separee du dispositif selon 1' invention. 

Connaissant le jeu de donn6es-image de I'objet 3D, et sa 
5 decomposition en mailles volumiques 14, ici t6tra6driques , le 
dispositif va determiner les forces dites "internes" c[ui sont 
exercSes sur chaque noeud 15 du maillage volumique. II est clair 
que la sorome des forces internes qui s'exercent sur un noeud 
donne n'est diff6rente de zero que lorsque ce noeud est mis hors 
10 de sa position d'6quilibre. Pour ce faire, il faut qu'une force 
externe soit appliqu6e sur I'un au moins des noeuds, par exemple 
par un outil, 

Les forces internes peuvent etre determin6es, soit sur la 
15 totality du maillage volumique & I'aide d'une premiere m^thode 
utilisant une unique technique dite "masses /tenseurs" qui sera 
decrite ci-aprfes, soit h I'aide d'au moins d'une seconde methode 
utilisant la premiere technique masses/tenseurs ainsi qu'une 
seconde technique, diff6rente> chaque technixjue 6tant appliquee 
20 k au moins un domaine issu d'une decomposition du maillage 
volumique. Get seconde methode, dite "hybride", sera decrite 
plus loin, en reference aux figures 14 et 15. 

Neanmoins, quelle que soit la methode, la determination des 
25 forces internes s'effectue grace & un module de forces internes 
16 (voir figure 2) qui est de preference integre dans la station 
graphique 2 (mais elle pourrait etre int6gr6e dans 1' unite de 
controle 1), et qui requiert 1' assistance d'un module de 
detection de collision 18 integre, de preference, dans 1' unite 
30 de controle 1. 

Un tel module a pour objet de determiner de fagon precise le 
lieu de collision (ou d' intersection) entre un point materiali- 
sant une partie de I'outil et I'enveloppe externe de 1' objet 3D, 
35 ou plus precisement une maille surfacique de cette enveloppe 
(voir figures 16 h 18). Cela peonnet de fournir au module de 
forces internes 16 la designation d'au moins un noeud de la 
maille surfacique sur lequel I'outil agit, et qui va permettre 
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de determiner les forces internes exerc^es en chaque autre noeud 
de la maille volumique. 

Le niveau de precision du point d ' intersection (on de collision) 
5 peut etre g6r6 h partir du nombre de points utilises pour 
mat^rialiser, d'une part, I'extr6mit6 7 de I'outil, et d' autre 
part, la tige 6 de celui-ci. En d'autres termes, la detection de 
collision effectu6e par le module de detection de collision 18 
peut etre appliqu6e en un ou plusieurs points mat^rialisant 

10 I'outil. Un exemple de discretisation de la forme de I'outil 
(tige et extr6mit6) est donn6 sur la figure 18, sur laquelle 
chaque point autre que les points mat6rialisant les noeuds d'une 
maille surfacique 13 sent relics les uns aux autres par des 
segments de fa<?on h reconstituer la forme de I'outil et plus 

15 particuli^rement celle de sa tige et de son extr^mite, ici des 
pinces coupantes . II est possible d'appliquer la detection de 
collision h la tige 6 de I'outil de fagon Si determiner une 
collision entre deux tiges d'outils actionn6s simultanement . 

20 Le module de detection de collision 18 fonctionne comme decrit 
ci-aprfes. Dans un premier temps, I'outil est discretise en un 
certain nombre de points pj^, qui depend 2l la fois de 1' instru- 
ment utilise et de la precision souhaitee- Puis, on d6finit un 
espace qui englobe I'objet 3D concern6. Get espace est alors 

25 subdivise en cellules volumiques, par exemple de formes 
paralieiepipediques . Le nombre de cellules volumiques dans 
chaque direction de 1' espace peut etre parametrable, et est de 
preference choisi en fonction de la geometric de I'objet 3D. De 
preference, le parametrage est choisi de sorte que le nombre de 

30 noeuds du maillage volximique de I'objet 3D soit sensiblement 
identique dans toutes les cellules volumiques de 1' espace defini 
par le module de detection de collision et contenant I'objet. On 
impose, 6galement de preference, une contrainte qui consiste h 
obliger que chaque cellule volumique intersectant la surface 

35 externe de I'objet englobe au moins un noeud du maillage. 

A partir de cette decomposition de 1' espace qui englobe I'objet 
3D en cellules vol\amiques, on g6nfere une table dans laquelle se 
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trouvent memorises des multiplets comportant les coordonn6es de 
chaque noeud du maillage volumique en reference ci la definition 
de la cellule volumique qui I'englobe. La determination d'un 
point de collision au niveau d'une maille surfacique 13 de 
5 I'enveloppe externe de I'objet 3D s'effectue k I'aide d'une 
comparaison entre la position d'un point mobile pj^ de I'outil 
qui est d6duite des informations fournies par I'interface- 
utilisateur h retour de force 4, ainsi qu'ci partir des multi- 
plets memorises dans la table. 

10 

Pour ce faire, on dSfinit une fonction de code qui permet de 
calculer une entree h l'int6rieur de la table et de trouver la 
liste des noeuds du maillage volumique qui correspondent ^ cette 
entr6e, lesdits noeuds appartenant alors k la meme cellule 
15 volumique que le point mobile pj^. Une telle fonction de code est 
donn^e par les formules referenc6es (1) en annexe, ci titre 
d ' exemple • 

On determine alors pour chaque noeud sii:u^ dans la meme 

20 cellule volumique que le point mobile py^, la distance eucli- 
dienne entre pj^ et chacun des Nj^. On note alors d^^^ la plus 
petite des distances ainsi calcul^es. On calcule ensuite la 
distance d qui s^pare le point mobile P]^ de chaque noeud Nj 
situ6 dans un nombre pr6d6teirmin6 de cellules volumiques 
25 voisines, par exemple 26. Si I'une de ces distances d associe k 
un noeud N est inf^rieure k dj^j^^' al<^rs d^^^^^ = d. II est clair 
que le temps de calcul sera d'autant plus important que le 
nombre de noeuds k l'int6rieur d'une cellule volumique sera 
grand • 

30 

Ensuite, on determine le nombre de mailles surfaciques 13 (ici 
des triangles) qui sont adjacentes au noeud N associ^ k la 
distance minimale d^,^ - d. Puis, pour chaque maille surfacique 
adjacente Si N, on cherche si elle est intersect^e par le segment 
35 [P^, ^t-il d^fini par les positions du point mobile pj^ aux 
instants t et t-1, comme illustrS sur la figure 17. Ainsi, 
quelle que soit la vitesse de d^placement de I'outil, on pourra 
tou jours determiner le point d' intersection (ou de collision). 
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6tant donn6 que seule la longueur du segment (P,., Pt-il varie en 
fonction de la vitesse. 

La position du point de collision est sub-triangulaire . Par 
5 consequent, les coordonn^es de ce point de collision sont 
d6tennin6es par une methode barycentrique . Ces coordonn^es 
barycentriques sont ensuite transmises au module de forces 
internes 16 de sorte que la force exerc6e par I'outil sur le 
maillage volumique soit estiinSe pr6cis6nient . 

0 

Connaissant la maille surfacique intersect6e, le module de 
forces internes 16 en d6duit les noeuds sur lesquels 1' action de 
I'outil s'exerce. Puis, a partir des coordonn^es barycentriques, 
il en d6duit la repartition du d^placement de I'outil virtuel 
15 sur chacun des noeuds de la maille surfacique d6sign6e. Dans le 
meme temps, 1' interf ace-utilisateur 4 fournit au module de 
forces internes 16 les informations relatives au d^placement 
vectoriel de I'outil, soit directement via un bus, soit via 
1 ' unit6 de controle 1 . 

20 

Le dispositif selon 1' invention utilise une technique que I'on 
appelle -masses/tenseurs" dans laquelle on soumet les noeuds 
appartenant h la maille surfacique, lieu de la collision avec 
I'outil, au d6placement vectoriel fourni (r^parti selon les 
25 coordonn6es barycentriques desdits noeuds), puis on estime les 
forces internes appliqu6es en chaque noeud du maillage volxami- 
que. Cette estimation est de preference iterative, chaque 
iteration permettant de calculer les forces internes appliquees 
en un premier "niveau" de noeuds, puis celles d'un "niveau" 
30 suivant constitue des noeuds voisins des noeuds de ce premier 
"niveau", et ainsi de suite jusqu'Si ce que tous les noeuds 
soient traites. Cela constitue une espfece de m6canisme de 
"propagation de contrainte" de noeuds St noeuds. En fait, lors de 
la premifere iteration seules les forces internes d'un premier 
35 niveau de noeud sont non nulles, celles des autres noeuds etant 
nulles. C'est au cours des iterations suivantes que progressive- 
ment ("niveau" aprSs -niveau") I'on obtient les forces internes 
(non nulles) des autres noeuds du maillage. 
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Cette technique permet d'estimer les forces internes de fagon 
ind6pendante de la direction des liaisons. Ainsi, la deformation 
induite par I'outil virtuel sur I'objet 3D est ind6pendante de 
sonmaillage volumique, Une telle caract6ristique d ' ind6pendance 
5 relativement au maillage initial, est particulierement int6res- 
sante dans le domaine de la simulation, notamment d' intervention 
chirurgicale, oil les deformations induites par un outil sont 
g6n6ralement volumique comme dans le cas d'une incision, ou d'un 
d^chirement, ou encore d'un retrait de matiere. 

10 

Dans une approche ou la loi de deformation est de type eiastique 
lin^aire volumique, la force interne exerc6e sur un noeud 
soumis ^ un d6placement u est donn6e par la formule r6f6renc6e 
(2) en annexe. 

15 

Cette formule ( 2 ) peut se decomposer en deux parties , ou deux 
forces, I'une associee au deplacement propre du noeud (ou masse) 
Ni par rapport h sa position d'6quilibre, 1' autre associ^e au 
deplacement d'un noeud voisin Nj. Soit v 1.' ensemble des noeuds 
20 Nj vols ins du noeud N^, la formule 2 des forces internes appl±- 
quees au noeud peut alors se reecrire comme indiqu6 dans la 
formule ( 3 ) donn^e en annexe . 

Dans cette formule (3), les expressions entre crochets, [T^^j^] et 
25 [Tij], representent des matrices de dimension 3x3 que I'on 
appellera par la suite respectivement tenseurs de noeud et 
tenseurs de liaison, du fait de leur definition similaire & 
celle des tenseurs de force bien connus de I'homme de I'art. 

30 Pour determiner ces tenseurs de noeud et de liaison, il est 
necessaire, tout d'abord, d'evaluer la matrice de rigidite 
eiementaire [K®] associee k chaque maille volumique (ici un 
tetraedre) adjacente & chaque liaison issue du noeud N^. 

35 Comme cela est represente sur la figure 7, le deplacement 
vectoriel d'un noeud entre un instant indice 0 et 1' instant 
t peut etre note Ax^^. Ce deplacement comprend trois composantes 
selon les trois directions de I'espace. 
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Une matrice de rigidity repr6sente pour un objet 3D, notaininent 
d^formable, les differentes relations qui existent entre ses 
616ments constituants . En d'autres termes, elle traduit 
1' influence d'un noeud sur d'autres noeuds, laquelle est donnee 
5 par la loi de deformation. Les valeurs num6riques des composan- 
tes d'une matrice de rigidite dependent par consequent des 
parametres d' elasticity de I'objet qui sont connus. La loi de 
deformation peut etre memorisee dans le module de forces 
internes 16 de la station graphique 2- 

10 

Si I'on designe par T une maille volumique tetr aedrique , alors 
la matrice de rigidite eiementaire de ce tetraedre est donn6e 
par la formule (4) donnee en annexe, & titre d'exemple. 

15 une fois la matrice de rigidite [K®] de chaque tetraedre T 
calcuiee, on associe ^ chaque sommet de tetraedre (ou chaque 
noeud), ainsi qu'^ chacune de ses liaisons, respectivement un 
tenseur de noeud T^^ ou un tenseur de liaison T^j, Qomme cela est 
illustre sur la figure 8. Cette association s'effectue ^ partir: 

20 des formules implicites donnees en annexe sous la reference (5). 

A partir de 1' ensemble des matrices de rigidite eiementaires 
[K^], on peut determiner la matrice de rigidite global [K] de 
I'objet 3D. une telle matrice est illustree, h titre d'exemple, 
25 sur la figure 9, dans le cas d'un objet 3D represente par 15 
noeuds references de 1 ^ F- Ici, les croix representent des 
valeurs non nulles, tandis que les points representent des 
valeurs nulles. 

30 Afin de simplifier les calculs des forces internes, il est 
possible de limiter le nombre de liaisons dans la maille 
volumique en filtrant les valeurs des composantes des matrices 
de rigidite eiementaires [K®], par exenple en les comparant h un 
seuil- Ainsi, toute valeur infer ieure h ce seuil est consider6e 

35 comme nulle, ce qui revient a supprimer des liaisons entre deux 
noeuds. Cette technique de simplification des matrices permet de 
ne determiner que les liaisons qui sont indispensables h la 
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transmission des contraintes (de deplacement. ) dans le maillage 
volumique. 

Par ailleurs, il est. 6galement possible de simplifier le calcul 
5 des tenseurs de liaison en tenant compte de la propri6t6 de 
sym^trie de la matrice de rigidity globale [K] . En effet, un 
seul et unique tenseur de liaison suffit pour caracteriser 
chaque liaison. En d'autres termes, soit l^j la liaison d6finie 
entre les noeuds et Nj^ le tenseur de liaison associe a cette 
10 liaison sera [T^j] si i>j et [T^j]^ si i< j . L'exposant t d6signe 
la transpos^e de la matrice concernSe. 

Arriv6 h ce stade, on dispose, d'une premiere part, des 
coordonn6es de chaque noeud N^^, d'une seconde part, de I'enve- 

15 loppe externe, du maillage surfacique de cette enveloppe externa 
et du maillage volumique final de 1 ' integrality de I'objet 3D, 
d'une troisi^me part des liaisons l^j entre deux noeuds Ni et Nj 
de 1' ensemble des t^traedres T composant le maillage volumique 
de 1' Ob jet 3Dr «t d'une qxiatri^me part,, de 1' ensemble des 

20 tenseurs de noeud [T^i] et tenseurs de liaison [T^j] associfis ^ 
chaque noeud et chaque liaison de chaque t6trafedre T. 

On d^finira alors, de pr6f6rence, une structure (ou base) de 
donnSes permettant de connaitre, d'une premifere part, I'indice 

25 des deux noeuds situ^s aux extr6mit6s de chaque liaison, d'une 
seconde part, la reference d'un t^trafedre (maille volumique) 
adjacent h un triangle (maille surfacique), d'une troisi^me 
part, l'6tat d'un t^trafedre (act if ou inactif), et d'une 
quatrifeme part, l'6tat d'une liaison (active ou inactive). Ces 

30 informations pourront etre utiles pour la prise en compte de 
deformations de type incision, d6chirement ou retrait de 
mati^re, aussi bien au niveau graphique que physique. Les 
diff^rentes Stapes d'obtention de tous ces paramdtres sont 
resumSes sur la figure 10. 

35 

Afin de tenir compte de forces de surface auxiliaires, permet- 
tant par exemple d' amplifier un effet ressenti lors d'une 
incision ou dScoupe, le module de forces internes 16 du 
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dispositif selon 1' invention peut combiner les forces de liaison 
surfaciques et volumiques decrites precedemment h des forces de 
tension y^^^^^^^^ du type de celles donnees dans la formule (6) 
en annexe. Ces forces peuvent etre choisies en faisant varier 1q 
5 la longueur & vide du ressort associ^, ainsi que sa raideur, 
selon 1' importance de 1 ' amplification souhaitee. Cette impor- 
tance peut etre definie, par exemple, par I'op^rateur, ou bien 
par I'interface-utilisateur 4 selon le type d'outil et/ou le 
type d'objet sur lequel 1 ' intervention a lieu. 

10 

Iiorsqu'une telle prise en compte est effectuSe, la formule (3) 
de la force interne exerc^e sur un noeud se transforme en la 
formule (7) donn^e en annexe. 

15 Par ailleurs, il est egalement possible de tenir compte dans le 
calcul des forces internes exerc6es sur chaque noeud du maillage 
volumique, d'un ou plusieurs types de forces ext^rieures dues h. 
la presence d'autres objets au voisinage de I'objet 3D concern^, 
ainsi qu'aux forces gravitationnelles <iui font qu'un objet 3D 

20 aura tendance h s'affaisser sous son propre poids du fait de la 
gravity. De telles forces ext6rieures peuvent Stre par exemple 
le poids exerc6 en une partie de 1' objet 3D par un objet voisin, 
ou bien la force d' attachement exercSe par un lien, comme par 
exemple un ligciment, entre 1 'objet 3D et un objet voisin. 

25 

II est Clair qu'une force de type attachement, peut etre 
consid6r6e coirane une condition aux limites assimilable, par 
exemple, St une force de rappel d'un ressort amorti poss6dant une 
extremity fixe, 1' autre extr6mit6 6tant relive 2i un noeud d'une 
30 maille surfacique de 1' objet 3D. 

La determination de la deformation d'un objet 3D, ou plus 
precis^ment de son maillage volumigue, requiert la connaissance 
chaque instant des d^placements de chaque noeud de ce maillage 
35 volumique. Cette determination en temps reel repose sur un 
principe ^ d' animation dynamique qui relie la somme des forces 
applic[uees sur 1 'objet 3D d6formable 2i son acceleration. Ici, 
1 'objet 3D deformable est forme par un ensemble de masses 
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ponctuelles (ou noeuds N^^) relives entre elles par des forces de 
liaison. II est par consequent possible d'appliquer 1' Equation 
g6n6rale de la m6canique a un maillage volumique de fagon 2l 
simuler le comportement dynamique de I'objet 3D qu'il represente 
5 lorsque celui-ci est d^form^ par des forces ext6rieures g^n6r6es 
par un outil dent le d6placement est d6fini par I'interface- 
utilisateur 2i retour de force 4, en reaction & une manoeuvre du ' 
joystick 5 par un op6rateur. 

10 A un instant t donne, la force appliqu6e sur un noeud de 
masse m^^ est done donn^e par 1' Equation g6n6rale de la dynamique 
donn^e dans la formule (8) en annexe. 

Dans la suite, on supposera que tous les noeuds du maillage 
15 volumique possddent la meme masse {m^ = ^/). Bien entendu, on 
pourrait proc6der dif f 6remment , chaque noeud poss6dant une masse 
propre, determinable par la connaissance de la constitution de 
la r6gion qu'il represente. Ces informations de constitution 
peuvent faire partie du jeu de donn6es-image formant 1' image 
20 num^rique 3D de I'objet. Mais elles peuvent Sgalement etre 
f oumies s^parSment . 

II est done possible de determiner 1' acceleration de chaque 
masse ou noeud N^, et par consequent de remonter dans un premier 

25 temps Si sa vitesse v^, puis dans un second temps a sa position 
^i+i ^ instant t+At en effectuant des integrations, par 

exemple num6riques, C'est ce calcul de toutes les positions x^^.^^ 
h 1' instant t-»-At qui va fournir les deplacements de chaque noeud 
et par consequent la deformation de la maille volumique. Pour ce 

30 faire, on pourra utiliser des methodes de type Euler, ou bien 
Runge-Kutta d'ordre 2, ou plus pr6ferentiellement d'ordre 4. 

La methode d' Euler est celle qui est la plus couramment 
employee. Elle consiste, par exemple, 2i partir d'une equation 
35 dif ferentielle dx / dt = f(x,t), avec comme condition initiale 
X = Xq & 1' instant t = 0. On choisit alors un pas temporal At 
tel que t^ = iAt, avec i = 0,1,2... Connaissant la position d'un 
noeud x. & 1' instant t., il est alors possible d'en deduire la 
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position de ce meme noeud ^ 1' instant, t^+j^ = t+At (voir formule 
( 9 ) en annexe ) . 

Deux additions^ deux multiplications et une Evaluation des 
5 forces sont alors n6cessaires pour obtenir la position d'un 
noeud N alnsi que sa vitesse ^ un instant t+At, a partir de son 
acceleration 1' instant t. La vitesse et la position du noeud 
sont donn6es & titre d'exemple dans les formules (10) en annexe. 

10 La formule de Runge-Kutta d'ordre 4 est actuellement pref^r^e en 
raison notamment de sa precision num^rique. Elle repose, par 
exemple, sur les formules (11) donnees en annexe. 

Dans cette m6thode, il est n6cessaire d'Evaluer quatre fois les 
15 forces externes et internes en chaque noeud. 

Bien entendu, d'autres techniques d' integration, notamment 
num^riques, peuvent Stre utilisees. 

20 Sur la partie gauche de la figure 11 se trouve repr6sent6 
1 ' enchainement des etapes qui permettent de determiner le 
deplacement de chaque noeud de la itiaille volumique, ou de la 
partie seiectionnee de celle-ci, h partir du deplacement 
vectoriel de I'extr^mite 7 de la tige 6 de I'outil virtuel 

25 relativement & I'objet 3D, du & 1' action de I'opErateur sur le 
joystick 5. Ce deplacement vectoriel est fourni par 1 'interface- 
utilisateur a retour de force 4, via le bus ISA, de preference 
en meme temps que la position de I'extremite 7. 

30 Des que le module de detection de collision 18 a determine un 
point d ' intersection entre "I'outil" et le maillage surfacique, 
le module de forces internes 16 peut commencer h determiner les 
forces internes exercees sur les noeuds de la premiere partie du 
maillage volumique. 

35 

Pour ce faire, il applique sur chaque noeud de la maille 
surfacique de collision detect6e (ici il y en trois du fait de 
la forme triangulaire de la maille) un deplacement vectoriel 
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egal ^ celui fourni par 1' interf ace-utilisateur 4 qui d6finit le 
d^placement de I'outil. Puis^ il determine la force interne 
exerc^e sur chaque noeud de la premiere partie du maillage en 
utilisant les formulas <3) ou (7) selon les forces externes 
5 choisies par I'operateur et/ou 1 ' interf ace-utilisateur 4. Cett-e 
determination fait appel aux tenseurs de noeud et de liaison, 
[T^il et [Tj^j], pr^c^demment determines et memorises. En fait, 
lors de la premiere iteration, bien que I'on calcule les forces 
internes appliquees en chaque noeud, seules celles des noeuds 
10 voisins des noeuds "de collision" sont non nulles, etant 

donne que les deplacements des autres noeuds sont nuls ci ce 
stade du calcul. On en deduit 1 ' acceleration , la vitesse et le 
deplacement des noeuds premiers voisins, et par consequent leurs 
nouvelles positions. 

15 

Ces nouvelles positions des noeuds deplaces sont alors figees 
pour la suite du calcul, et le module de forces internes 16 
reitere son calcul de fagon ei determiner en chaque noeud de la 
premiere partie <ici toute la -ma^ille) , la force interne (non 
20 nulle) qui lui est appliquee, et par consequent son deplacement. 
Le calcul des forces internes s'est done effectue selon un 
modeie de type "propagation de contraintes" • On connalt alors la 
deformation induite par le deplacement vector iel de "I'outil" 
virtuel sur la maille volximique representant I'objet 3D. 

25 

Entre deux boucles de calcul des forces internes et des 
deplacements dus h un deplacement vectoriel de "I'outil", c'est 
k dire en cours de propagation de contraintes^ il est particu- 
lierement avantageux de tenir compte d' informations transmises 

30 par 1' interf ace-utilisateur 4, ou definies par I'operateur, 
portant sur le type d'outil manoeuvre et par consequent le type 
d' intervention effectuee. En effet, de telles informations 
peuvent indiquer au module de forces internes 16 que 1' action de 
I'outil est destinee & couper, ou inciser, ou supprimer de la 

35 matiere. Ainsi, apres une iteration au cours de laquelle on a 
determine des forces internes (non nulles) exercees sur certains 
noeuds, le module de forces internes 16 peut supprimer les 
liaisons qui sont dans le plan de coupe d'un outil d' incision ou 
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retirer la matifere (mat6rialis6e par des mailles volumiques, ici 
des t6tra6dres) qui se trouve dans la zone d' interaction de 
"I'outil" virtual, sur la base des informations d' action qui lui 
sont fournies par le module de detection de collision 18 et 
5 1' inter face-utilisateur 4. 

Dans ce cas, un certain nombre de liaisons et/ou de noeuds sont 
supprim6s. Le module de forces internes doit alors mettre ^ jour 
les tenseurs de noeud et les tenseurs de liaison en fonction de 
10 la ou des liaison(s) et/ou noeud(s) supprim^s. Cette prise en 
compte des informations pr6cit6es tient lieu de critdre de 
suppres s ion . 

Les outils virtuels capables de d6coupe (ou d' incision) et/ou de 
15 retrait de matifere sont par exemple des scalpels, pinces 
coupantes, ou bien bistouri m6canique ou 61ectrique, ou bien 
encore lasers. 

Par ailleurs, il est 6galement particuli^rement ava:ntageux 
20 qu'entre deux boucles de calcul le critfere de suppression soit 
capable d'effectuer un test portant sur I'int6grit6 des liaisons 
entre noeuds dont les dSplacements viennent d'etre d6termin6s et 
sont d^sormais fig^s. 



35 



25 ce test a pour but de determiner si certains d6placements sont 
tels que dans la r6alit6 ils auraient caus6 la rupture d'une ou 
plusieurs liaisons et/ou la suppression d'un ou plusieurs noeuds 
et done d'un ou plusieurs t6tra%dres (ou mailles volumiques). Le 
test consiste en la comparaison des d6placements h un seuil de 

30 rupture de liaison pr6d6termin6 . Ce test pourra porter, par 
exemple, sur la variation de volume de la maille volumique 
comprenant la maille h supprimer, et/ou la variation de longueur 
des liaisons de la maille volumique comprenant la liaison a 
supprimer . 



Si aucune liaison ne doit §tre supprim^e, le module de forces 
internes 16 passe Si 1' iteration suivante (boucle suivante) de 
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fagon k poursuivre le calcul des forces internes (poursuite de 
la propagation de contraintes ) . 

En revanche^ si le test indique qu'au itioins une liaison doit 
5 etre detruite, c'est-^-dire si 1' intensity du d^placement 
dSpasse le seuil predetermine, alors le module de forces 
internes 16 proc^de k une mise k jour des tenseurs de noeuds 
[T^j^] et des tenseurs de liaison [T^^j] de la maille volumique ou 
de la premiere partie de celle-ci, en tenant compte des 

10 d^placements calculus. La valeur de certains tenseurs de noeud 
corome de liaison peut alors devenir nulle, C'est notamment le 
cas, lorsque les quatre liaisons qui relient un noeud au quatre 
autres noeuds d'un t6tra6dre donne se trouvent rompues. Dans ce 
caSr il faut supprimer un noeud de la maille volumique, ou 

15 plusieurs si cela s' impose, ainsi que toutes les liaisons 
associees . 

L' operation ne s'arrete pas a ce stade, 6tant donn^ que, pour 
des raisons d'fiquilibre de la" maille "volttmiq^^ il faut 

20 que cette dernidre conserve sa structure g^n^rale de depart au 
niveau de ses maillages surfacique et volumique. Or, le retrait 
d'un noeud, ou la rupture d'une liaison provoquent des modifica- 
tions d'une ou plusieurs mailles, celles-ci n'^tant plus, dans 
1' example illustre, tetraSdriques ou triangulaires selon que 

25 I'on parle de maille volumique ou de maille surfacique. Un tel 
cas est illustre sur la figure 12 sur laquelle 1' action de 
I'extr6mit6 7 de I'outil a provoqu6 la rupture de quatre 
liaisons suivant un plan de coupe particulier. Chacune des 
quatre liaisons se trouvant divis^e en deux, il est alors 

30 n6cessaire de rajouter des noeuds, et par consequent des 
liaisons, de fagon h, reconstituer des mailles t6tra6driques dans 
la zone incis6e. En d' autres termes il faut remailler localement 
le maillage volumique. 

35 £tant donne que I'ajout de noeuds a pour consequence d'accroltre 
la masse totale de I'objet 3D, il est alors necessaire de 
modifier la masse de tous les autres noeuds du maillage 
volumique, ce qui peut etre fait par une repartition homogene. 
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De pluSr I'ajout de noeuds entraine I'ajout de liaisons pour que 
les mailles conservent leurs formes g^n^rales, ici triangulaires 
ou t§tra6driques • II faut. par consequent que le module de forces 
internes 16 cr^e et/ou d^truise des matrices de rigidity 
elementaires [K®] et des tenseurs de noeud et de liaison^ et 
recalcule les matrices de rigidity 616mentaires et les tenseurs 
de noeud et les tenseurs de liaison pr6c6demment calculus, de 
sorte que la poursuite des calculs des forces internes exercSes 
sur les noeuds puisse se faire en tenant compte de ces modifica- 
tions • 

II est Clair qu'il est particulierement avantageux que le module 
de collision 18 puisse mettre a jour ses multiplets au fur et a 
mesure de 1' Evolution (deformation) de la partie de la maille 
volumique . 

D6s que le calcul des forces internes et d6placements est 
consider^ par le module de forces internes corame termini pour la 
valeTir de d^placement vectoriel de I' outil fournie, le module de 
rafraichissement 17 procSde ^ la determination des nouvelles 
donnees-image de I'objet et de la materialisation de I'outil en 
vue de leur affichage^ tandis que dans le meme temps, le module 
de forces internes recommence un nouveau calcul avec un nouveau 
deplacement vectoriel de I'outil appliqu6 eventuellement a 
d'autres noeuds que ceux sur lesquels ils ont et6 initialement 
appliques • 

Sensiblement simultanement h la determination des nouvelles 
donnees-image de I'objet et de "I'outil-, un module de reaction 
20, de preference integre dans 1 'unite de controle 1, determine 
la force de reaction de I'objet, laquelle correspond 2i sa 
deformation (deplacements des noeuds) estimee h partir des 
forces internes. Pour chaque noeud d6place par I'outil, on 
connait sa position initiale (St t=0) et sa position en cours (& 
1' instant t) - On connait par ailleurs les tenseurs associes h ce 
noeud et les liaisons qui lui sont associees, ce qui permet 
d'appliquer la formule (3) donnee en annexe. 
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La force de reaction est alors transmise h 1 ' interf ace-utilisa- 
teur ^ retour de force (LIE) de sorte qu'elle g6n6re sur le 
joystick force (ou retour de force) sensiblement 6quilibr6e par 
ladite force de reaction. Cela pentiet de transmettre & l'op6ra- 
5 teur des reactions physiques de I'objet, ce qui est essentiel ^ 
un bon controle de la manoeuvre d' intervention, 

Coirane il ressort de la description ci-avant, la technique de 
calcul des masses/tenseurs permet de simuler tous types de 
10 deformations^ qu'elles soient de type g6om6trique, de type 
incision ou d^coupe, de type d6chirement ou fracture, ou bien 
encore de type retrait de inatiSre, 

On a repr6sent6 sur la figure 13 quatre images d'un foie humain 
15 soumis ^ 1' action d'un outil virtuel de type pinces coupantes, 
Ici 1' action est une simple compression. 

Afin d'am61iorer la quality de 1' image, on peut envisager, avant 
chaque mise h jour de la sc^ne af f ich6e> <ie r^it^rer un c?ertain 

20 nombre de fois le calcul complet de la deformation de la maille 
volumique, et par consequent les tests portant sur les suppres- 
sions de noeud et de liaison, ainsi que ceux portant sur les 
valeurs des composantes des tenseurs. Ainsi, tant que le nombre 
d' iterations est infer ieur h un seuil choisi, et que le temps de 

25 calcul reste inferieur au temps necessaire au raf raichissement 
des Images sur les moyens d'affichage, alors on peut recommencer 
la boucle* En revanche, si tel n'est pas le cas, on procfede ^ la 
mise h jour de la scene affichee. 

30 Toutes les principales 6tapes de la technique masses/tenseurs 
sont illustrees dans le diagramme de la figure 11* 

L' invention propose egalement un precede eiectronique de 
traitement de donnees-image pour la mise en oeuvre du dlspositif 
35 precite, dans lequel on prevoit les principales etapes suivan- 
tes: 
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- tout: d'abord, pr6voir une interf ace-utilisateur 4 capable de 
donner un retour de force ^ conf orm^ment aux reactions d'un 
outil t 

- puis, estimer un point d' intersection (ou de collision) entre 
5 une droite materialisant un deplacement tir6 de 1' action d6finie 

par 1 ' interf ace-utilisateur et le itiaillage surfacique, 

- puis, 6tablir un champ de forces internes entre des noeuds 
d'une premifere partie au moins d'un maillage volumique fonction 
d'un maillage surfacigue d'un objet 3D apparaissant dans un jeu 

10 de donn6es- image, & partir d'une loi de deformation , par exemple 
de type ^lastique linSaire voluiaigue, d'un deplacement, induit 
par une action d^finie par 1 ' interf ace-utilisateur et repre- 
sentative d'une manoeuvre de 1' outil, applique aux noeuds 
appartenant la maille surfacique contenant le point d'inter- 

15 section, de conditions aux liiaites, et de tenseurs de noeud et 
de tenseurs de liaison issus respectivement pour chaque noeud et 
chaque liaison de cette premiere partie au moins, de matrices de 
rigidite propres ci chaque maille volumique d'au moins la 
premidre partie et fonction de La Xoi de deformation, 

20 - ensuite, determiner la force de reaction de 1' objet qui 
correspond & sa deformation estimee & partir des forces 
internes, de sorte que la force g6n6r6e par 1' interf ace-utilisa- 
teur soit sensiblement equilibree par cette force de reaction, 

- et enfin, calculer de nouvelles donnees-image de 1' objet, en 
25 presence des deformations estimees, compietees de la representa- 
tion de 1' outil. 

Dans certains cas complexes, ou le nombre de donnees-image 
formant 1' image numerique tridimensionnelle d'un objet 3D est 

30 trfes important, le temps necessaire au calcul des forces 
internes et de la deformation de 1' integralite du maillage 
volxomique peut devenir superieur au temps de raf raichissement 
des images sur les moyens d'affichage (ici un moniteur). II en 
resulte que 1 'application de la technique decrite ci-avant ne 

35 peut etre faite sur 1 ' integralite de 1' objet 3D. C'est pourquoi 
le dispositif selon 1' invention peut appliquer aux donnees-ima- 
ges une seconde methode dite hybride dans laquelle la technique 
des masses /tenseurs ne va etre appliquee qu'& une preiciiere 
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partie du maillage volumique, tandis que la partie compl6mentai- 
re, dite seconde partie, de ce maillage volumique va etre 
trait^e ci I'aide d'une seconde -technique de calcul en temps r6el 
reposant par exemple sur une approche de type ^ISments finis, 
faisant appel h des pr6-calculs m6moris6s. Cette seconde 
technique est notamment d^crite dans 1' article d4j^ cit6 de S. 
Cot in, H. Delingette, M. Bro-Nielsen et N. Ayache, "Geometric 
and physical representations for a simulator of hepatic 
surgery" , publi6 dans les acties de la Conference Medecine meets 
with virtual reality de janvier 1996. 

Dans ce qui suit, on appellera premiere partie la partie du 
maillage volumique ^ laquelle est appliqu6e la technique des 
masses /tenseurs , et seconde partie la partie compl6mentaire de 
la premiere partie, h, laquelle est appliqu6e la technique des 
616ments finis via des pre-calculs m6moris6s . 

Dans le module hybride on subdivise le maillage volumique 
(initial ou final) en au moina <ieux <lomaines. II est clair que 
I'on pourrait envisager une subdivision en trois domaines, voire 
plus, grlice & un traitement, par exemple en paralldle, des 
diff brents domaines. Cette subdivision en domaines peut etre 
effectu4e par le dispositif selon 1' invention, lorsque celui-ci 
comporte un module de maillage sp6cialis6 du type de celui 
d^crit pr6c6demment , mais il peut 6galeitient etre effectu6 en 
dehors du dispositif dans le cas contraire. 

Qu'elle soit faite par un module externe ou bien par un module 
interne au dispositif, la subdivision en domaines r6sulte d'une 
selection par un op^rateur ou d'une detection automaticjue 
resultant d'un critfere predetermine applique aux donnees-image 
de I'objet 3D. Le critere pourra etre notamment une comparaison 
d'intensite entre les differentes donnees-image ou bien entre 
ces donnees-image et un seuil, par exemple dans le cas d'une 
tumeur, ou une detection d'un parametre physique ou anatomique 
contenu dans lesdites donnees-image de I'objet 3D, et designant 
des zones complexes dans lesquelles se trouvent des sous- 
structures corome par exemple des tumeurs. 
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On peut envisager que le maillage soit subdivis6 en une 
multiplicity de r6gions, par exemple six ou sept, et que celle 
dans laquelle est d^tectee la collision avec "I'outil" soit 
d6sign6e automat iquement comma premiere partie, la partie 
5 compl6mentaire formant alors la seconde partie du maillage 
volumique . 

Le dispositif pourra comprendre un module sp6cifique de 
partitionnement, different du module de maillage, pour d6termi- 
10 ner la subdivision en domaines. Celui-ci sera de pr6f6rence 
int6gr6 dans le module de forces internes 16. 

Dans le module hybride, une fois la subdivision en parties 
effectu6e, le module de forces internes 16 debute le calcul des 

15 forces internes de 1 ' integrality du maillage volumique de 
I'objet 3D. Ce calcul peut d^buter soit par la premiere partie 
du maillage, & I'aide de la technique masses/tenseurs, soit par 
la seconde partie du maillage, k I'aide de la technique des 
yiyments finis via des pr6-ealcnils (technique dite ^des pr6- 

20 calculs") qui repose sur deux principes dits de superposition et 
de linearity. 

Le d^placement vectoriel de I'outil induit un d6placement 
vector iel de certains noeuds du maillage volumique, lequel est 

25 considers comme un multiple d'un deplacement vectoriel yiemen- 
taire. Puis, le d6placement total de chaque autre noeud i (ceux 
sur lesquels Inaction de I'outil n'est pas directement exerc^e) 
est d6tenain6e en effectuant une sommation (principe de 
superposition) des contributions des d^placements de tous les 

30 noeuds j ayant subis un deplacement. 

Afin de mettre en oeuvre ce principe de superposition, il faut 
determiner une table de pr6-calcul de n x n tenseurs de 
deformation, oCi n represente le nombre total de noeuds du 
35 maillage, ou d'une partie de celui-ci (par exeit^jle uniquement 
les noeuds de la surface). Pour ce faire, on applique un 
deplacement eiementaire en un noeud donne, puis on calcul la 
deformation que ce deplacement induit sur tous les autres noeuds 
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du maillage. On en d6duit n premiers tenseurs de d6forination 
6 lament aires , un pour chaque noeud- Puis on recommence les deux 
operations pr§cit6es en appliquant le d^placement 616ment.aire k 
un autre noeud du maillage, ce qui fournit n seconds tenseurs de 
5 deformation. On renouvelle ces deux operations pour les n noeuds 
du maillage, jusqu'Si obtention des n x n tenseurs de deformation 
que I'on memorise sous forme de table dite de pre-calcul. 

On dispose alors d'une source de renseignements qui permet de 
10 savoir instantanement la deformation induite sur la maille par 
un deplacement eiementaire d'un noeud donn6. II n'y a plus alors 
qxi'h determiner le rapport entre le d6placement vectoriel du 
noeud et le deplacement eiementaire, puis St multiplier le 
deplacement des noeuds du maillage memorise dans la table par ce 
15 rapport (principe de linear ite). En d'autres termes, le 
deplacement total d'un noeud est obtenu en sommant les 

contributions de tous les tenseurs de deformations T'^j designes 
par I'indice i dans la table (T ' T' ^^^T' • • • ' in) sans 
oublier de multiplier ces contributions par le rapport. 

20 

Des elements descriptifs plus detailies de la seconde technique 
dite des pre-calculs sont indiqu6s, notamment, dans la publica- 
tion de S.Cotin, H-Delingette et N.Ayache, "Real-time non linear 
elastic deformations of soft tissues for surgery simulation" , 
25 Rapport de recherches de I'INRIA, qui sera rendue publique apr&s 
dep6t de la presente Demande de Brevet. 

Bien entendu, comme dans la technique des masses /tenseurs, la 
seconde technique tient compte des conditions aux limites, et 
30 autres forces auxiliaires et/ou externes. 

Quelle que soit la partie du maillage par laquelle d6bute le 
calcul des forces internes, 11 faut tenir compte des noeuds "de 
connexion- 19 (voir figure 14) qui se trouvent plac6s Si I'inter- 
35 face entre les deux parties, et qui par consequent appartiennent 
en meme temps ^ ces deux parties. Selon 1' invention, les noeuds 
de connexion d6finissent des conditions aux limites et par 
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consequent ils sont des parametres de la deformation du maillage 
complet • 

Si I'on raisonne a partir d'un debut de determination des forces 
5 internes sur la premiere partie I'aide de la technique des 
nvasses/tenseurs ) ^ mais on pourrait faire I'inverse, on determine 
les forces internes^ et par consequent les d^placements, de tous 
les noeuds de la premiere partie du maillage, y compris ceux des 
noeuds de connexion. Les d6placements de ces noeuds vont alors 

10 permettre de determiner h I'aide de la table n x n de tenseurs 
de deformation T'^j pre-calcuies et memorises, selon la seconde 
technique decrite precedemment , les deplacements de tous les 
noeuds de la seconde partie du maillage • Cela fournit de 
nouvelles conditions aux limites (forces externes) pour la 

15 premiere partie du maillage, lesquelles vont etre prises en 
compte pour recalculer les forces internes appliquees en chaque 
noeud de la premiere partie, et par consequent les deplacements 
de ses noeuds, y compris ceux des noeuds de connexion, ce qui 
foumiti de nouvelles conditions aux limites pour la 4&econde 

20 partie du maillage qui vont permettre h leur tour de recalculer 
les deplacements des noeuds de la seconde partie • 

On reitere cette procedure, de preference, jusqu'k obtention 
d'un equilibre entre les deux parties du maillage, au niveau des 

25 forces exercees sur les noeuds de connexion. Get equilibre peut 
etre fixe par un seuil de convergence, ou bien par un nombre 
d' iteration "flottant" dependant du temps de calcul autorise, 
lequel doit etre imperativement inferieur au temps necessaire au 
raf raichissement d'une image (frequence superieure k 20 Hz). La 

30 boucle d' animation du modeie hybride est illustree sur la figure 
15 ou les maillages Ml et M2 representent respect ivement les 
premiere et seconde parties du maillage volumique. 

II est Clair que le module de forces internes 16 tient compte 
35 des suppressions et ajouts de liaisons et/ou de noeuds lorsqu'il 
utilise le modSle hybride. Par consequent , les tenseurs de noeud 
et de liaison utilises dans la premiere technique masses/ten- 
seurs sont reactualises en temps reel. 
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D^s que le calcul des forces internes et d6placeinents est 
consid6r6 par le module de forces internes 16 coirane termini pour 
une valeur de d^placement vector iel de I'outil^ ledit module de 
raf ralchissement 17 proc^de ^ la determination des nouvelles 
5 donn^es-image de I'objet et de la materialisation de I'outil en 
vue de leur affichage, tandis que dans le meme temps, le module 
de forces internes 16 recommence un nouveau calcul avec un 
nouveau d6placement vectorial de I'outil appliqu6 6ventuellement 
Si d'autres noeuds que ceux sur lesquels ils ont 6t6 initialement 
10 appliques. 

II est Clair que dans le cas d'un modele hybride, il est 
particuliferement avantageux que le module de collision 18 
d6finisse des sous-tables de multiplets pour chaque partie de la 
15 maille volumique. Par ailleurs, il est particuliferement 
avantageux, dans un souci de precision dans la detection des 
points de collision, que les multiplets des sous-tables soient 
mis h jour au fur et a mesure de 1' evolution (deformation) des 
differentes parties de la maille volumique. 

20 

Le dispositif selon 1' invention peut etre implante dans une 
memoire, par exemple de masse, d'un ou plusieurs moyens de 
calcul de type station de travail et/ ou ordinateur, sous, la 
forme de modules logiciels. 

25 

A toutes fins utiles, il est precise que des elements descrip- 
tifs plus detailies sont indiques dans la these de doctorat de 
S.Cotin, "Modeies anatomiques deformables en temps-reel - 
Application h la simulation de chirurgie avec retour d' effort", 
30 soutenue le 19 novembre 1997, ainsi que dans la publication d6ja 
citee "Real-time non linear elastic deformations of soft tissues 
for surgery simulation", de S.Cotin, H.Delingette et N.Ayache. 
Ces documents seront rendus public apres depot de la presents 
Demande de Brevet. 

35 

L' invention ne se limite pas aux modes de realisation deer its 
ci-avant, seulement h titre d' exemple, mais elle englobe toutes 
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les variantes que pourra envisager I'homme de I'art dans le 
cadre des revendications ci-apr6s. 

Ainsiy on a d6crit. un dispositif de traitement d' image 3D, et le 
5 proc^dS associe, dans lesquels 1' action sur un outil 6tait 
induite par la manoeuvre d'un "joystick" par un op^rateur, puis 
d6finie par 1 ' inter face-utilisateur. Mais, on peut envisager que 
1' action sur 1' outil soit directement d6finie par I'interface- 
utilisateur, cette action 6tant mSmoris^e et decomposable en 
10 sous-actions 616mentaires , l'op6rateur n ' intervenant plus que 
pour decider du choix de certaines sous-actions 616mentaires, 
par exemple en vue d'une formation (ou apprentissage) • 

Par ailleurs, on a d^crit des dispositif s et le proc6d6 associ^ 
15 destines a traiter des images en particulier m6dicales, et 
notamment du foie, Mais il est clair que 1' invention n'est pas 
limit6e k ce seul domaine. Elle s ' applique ^galement a d'autres 
secteurs du domaine medical, ainsi qu'aux domaines non mSdicaux 
pour lesquels le traitement en t^mps r^l d' images -d'objets 3D 
20 sur lesquels une simulation d' intervention, au sens le plus 
large du terme, revSt un int6r6t particulier aussi bien sur le 
plan de I'enseignement des techniques d ' intervention que sur le 
plan du perfect ionnement ou de la mise au point de nouvelles 
techniques d ' intervention . 

25 

De plus on a d^crit des modules de traitement dans lesquels les 
matrices et tenseurs 6taient calculus ci partir d'une loi de 
deformation 61astique lin^aire volumique. Cependant, d'autres 
types de lois de deformation peuvent etre utilises, et notamment 
30 des lois de type non lineaire, 

Enfin, on a d6crit un module "hybride" faisant appel deux 
techniques diff ^rentes dont une premiere technique de calcul en 
temps reel appeiee "masses /tenseurs ** , la seconde technique etant 
35 dans I'exeii^le choisi une technique de calcul en temps reel dite 
"des pre-calculs" reposant sur une approche par elements finis. 
Mais, il est clair que toute autre seconde technique pourra etre 
envisagee, qu'elle fasse appel ou non k des pre-calculs, pourvu 
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qu'elle permetlie un traitement des donnees-image compatible avec 
une simulation dynamique d'une intervention en temps r6el. 
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Annexe 

* Formules ( 1 ) : 
5 code = X + y * N^eiis x ^ * ^cells y * ^cells x 



, Px - bmin-. 



^^■'■■'■wielth y 



. ^, PZ - bZttin^ ^ 



- "int" : fonction donnant la partie entifere de son argument, 

- biain^ <i=x,y,z) : bornes inf^r inures de I'espace englobant 
10 I'objet, 

- cell^i^th i (i=x,y,2) : ddLmensions d'une cellule volumique, 

- N^eiis i (i=X'y) • nombres de cellules volumiques dans les 
directions de I'espace x et. y. 

15 * Formule (2\ : 

- u : vecteur de dimension 3*(N^+1) repr6sentant le d^placement 
20 du noeud et de ses vols ins, 

- : nombre de voisins du noeud N^, et 

- [M] : matrice de d^foirmation de dimension 3*3*(N^+1) repr6- 
sentant la loi de deformation de I'objet (connue par ailleurs). 

25 * Formule (3) : 
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* Formule ( 4 ^ 
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[K^ = Volume (T) [r^]nE^] IT^] 

(Ey] est, par exemple, donnS par la matrice : 



10 



[Et] = 



15 



A+2^i 


0 


0 


0 


A 


0 


0 


0 


A 


0 


U 


0 




0 


0 


0 


0 


0 


0 


0 


u 


0 


0 


0 


fj 


0 


0 


0 




0 




0 


0 


0 


0 


0 


A 


0 


0 


0 


A+2^ 


0 


0 


0 


A 


0 


0 


0 


0 


0 




0 


u 


0 


0. 


0 




0 


0 


0 


^ 


0 


0 


0 


0 


0 


0 


0 




0 


li 


0 


A 


0 


0 


0 


A 


0 


0 


0 


k+2fj 



20 



A et ^ : coefficients de LAM£, 

[r.p]^ est donnS, par example, par la matrice 



25 



[r„]t= 



30 



35 



ail 


al2 


a22 


0 


0 


0 


0 


0 


0 


a21 


a22 


ol23 


0 


0 


0 


0 


0 


0 


a31 


a32 


a33 


0 


0 


0 


0 


0 


0 


a41 


a42 


a43 


0 


0 


0 


0 


0 


0 


0 


0 


0 


ail 


al2 


a22 


0 


0 


0 


0 


0 


0 


a21 


a22 


a23 


0 


0 


0 


0 


0 


0 


a31 


a32 


a33 


0 


0 


0 


0 


0 


0 


a41 


a42 


a43 


0 


0 


0 


0 


0 


0 


0 


0 


0 


all 


al2 


a22 


0 


0 


0 


0 


0 


0 


a21 


a22 


a23 


0 


0 


0 


0 


0 


0 


a31 


a32 


a33 


0 


0 


0 


0 


0 


0 


a41 


a42 


a43 



- Ojj^ : premieres dSriv^es partielles des coordonn^s barycentri- 
ques Aj(x) d'un point x situ6 k l'int6rieur d'un tdtra^dre T : 
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^^^^^^ = o.ii, i = 1,2,3; j = 1,2,3,4 
OX,- 



* Formules ( 5 \ : 



5 a) pour le tenseur de force [T^j^] : 



Tfc^j + 52 "^^V^ic, s^*4i ; ic - 1,2,3; I =1,2,3 



- = local (m^) : indice local du noeud N^r 
10 - "local" : fonction de transfer!: affect ant un indice local 
(compris entre 1 et 4 dans le cas d'un t6tra6dre et correspon- 
dant a I'un des quatre sommets de T) au noeud represents ici 
par sa masse m^^. 

15 b) pour le tenseur de liaison [T^^j] : 



^Jt,i*^ ^Jc,J 53 ^Si^^k.sj*^! /' -fc"" 1#2,3; 1 - 1,2,3 



- s^ = local (la^) et Sj = local (nij) 



20 * Formula ( 6 ) : 



{■1.0 J) «V 



Fi =0 sxnon 



- < 1q, ou Iq est la longueur h vide d'un ressort de raideur 
choisie selon 1' importance de 1' amplification souhaitSe. 



BNSDOCID: <FR___2771202A1J_> 



2771202 

44 



* Formule ( 7 \ : 



5 * Formule ( 8 ) : 
* Formule ( 9 ) : 



10 



15 



20 



25 



Xi+1 = Xi + At*i'(x_i, ti) 

* Formules flO^ : 
v^.i = + At * 

= Pi + At * 

- : posit.ion du noeud N i 1' instant t^. 

* Formules (11) : 



fixj = At.f (x^, tp 

SXj = At.f (Xi + ^6xi, t^i-^At) 

&X3 = At.f (Xi+-|6xi, ti+-|At) 

6x4 = At.f (Xi+SXj, t^+AO 

^i.i = x,-.|&x,*|&x,*-|6x3*l6x. 



30 - la Vitesse v(t + At) = Vi^i at la position P(t + At) - Pi+i se 
d6duisent de 1 ' acceleration ai par les Equations suivantes : 



Svj^ = At. 
= At. Vi 



35 
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6-2 


= At. 






= At. 


" 1 




= At. 




6P3 


= At. 


V'i 




= At. 


a'" 


5P, 


= At. 


v"' 



- a^' , a^'' et a^''' : accelerations calcul^es respectivement 
partir des positions "interm^diaire" 5?^, 5P2 et 6P3. 
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1. Dispositif 61ect:ronique de traitement de donn6es- image, du 
type comprenant : 
5 - une interf ace-utilisateur (4) capable de g6n6rer un retour de 
force, conform6inent aux reactions d'un outil (6,7), 

- un module de forces internes (16) apte, sur designation d'un 
objet 3D apparaissant dans un jeu de donn6es-image, li 6tablir un 
champ de forces internes entre des noeuds (N) d'un maillage 

10 volximique fonction d'un maillage surfacique de cet objet, ^ 
partir d'une loi de deformation et d'une action d^finie par 
I'interface-utilisateur (4) et representative d'une manoeuvre 
dudit outil, 

- un module de reaction (20), pour determiner la force de 
15 reaction de 1' objet qui correspond a sa deformation estimee h 

partir des forces internes, de sorte que la force generee par 
I'interf ace-utilisateur (4) soit sensiblement equilibree par 
cette force de reaction, 

- un module de rafraichissement <17)^ pour calculer de nouvelles 
20 donnees-image de 1 'objet, en presence de ses deformations 

estimees, compietees de la representation dudit outil, 

caracterise en ce qu'il comprend un module de "collision" 
(18) propre h estimer un point d' intersection entre une droite 
materialisant un deplacement tire de ladite action d6finie et 

25 ledit maillage surfacique, et 

en ce que le module de forces internes (16) est agence pour 
estimer la force interne exercee sur chaque noeud d'une premiere 
partie au moins du maillage volumique de 1 'objet ^ partir dudit 
deplacement applique aux noeuds appartenant St la maille 

30 surfacique contenant ledit point d' intersection, de conditions 
aux limites, et de tenseurs de noeud et de tenseurs de liaison 
issus respectivement pour chaque noeud et chaque liaison de 
cette premiere partie au moins, de matrices de rigidite propres 
k chaque maille volumique d'au moins ladite premiere partie et 

35 fonction de la loi de deformation. 

2. Dispositif selon la revendication 1, caracterise en ce 
qu'il comprend un module de maillage () propre h designer ledit 
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Ob jet 3D par une determination d'enveloppe externa^ puis h 
decomposer ladite enveloppe en lesdites mailles surfaciques en 
vue de la decomposition du volume interne de 1' enveloppe en 
lesdites mailles volximiques. 

5 

3. Dispositif selon I'une des revendications 1 et 2, caract6- 
ris6 en ce que ladite enveloppe externe est obtenue par une 
methode de segmentation, en particulier par extraction d'iso- 
surf aces • 

10 

4. Dispositif selon I'une des revendications 2 et 3, caract6- 
rise en ce que lesdites mailles surfaciques sont des triangles. 

5. Dispositif selon la revendication 4, caract6ris6 en ce que 
15 lesdites mailles volumiques sont des tetrafedres formes h partir 

des mailles surfaciques triangulaires . 

6. Dispositif selon la revendication 5, caract6ris6 en ce que 
lesdlt.es mailles volumiques sont obtenues par une m^thode de 

20 type Delaunay-Voronoi. 

7. Dispositif selon I'une des revendications 2 St 6, caract6- 
ris6 en ce que le module de forces internes (16) est agenc6 poxir 
calculer les matrices de rigidity de chaque maille volumique, 

25 ainsi que lesdits tenseurs de noeud et tenseurs de liaison. 

8. Dispositif selon I'une des revendications 1^7, caract6- 
ris6 en ce que ladite loi de deformation est de type eiastique 
lineaire volumique. 

30 

9. Dispositif selon I'une des revendications 1 a 8^ caract6- 
ris6 en ce que le module de forces internes (16) est agenc6 pour 
calculer la force interne exerc6e sur chaque noeud de ladite 
premiere partie du maillage volumique & partir du produit de son 

35 tenseur de noeud et du d6placement estim6 de ce noeud, et d'une 
sommation, sur 1' ensemble des noeuds voisins poss6dant une 
liaison avec ledit noeud, du produit du tenseur de liaison. 
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associ6 h la liaison entre le noeud voisin et ledit noeud, et du 
d6placement estiin6 de ce noeud voisin. 

10. Dispositif selon I'une des revendications 1 Si 9, caract6- 
ris6 en ce que ledit module de forces internes (16) est propre 
h determiner les forces internes exerc^es sur certains au moins 
des noeuds de la premiere partie du maillage volumique & partir 
de la loi de deformation stockSe et de forces auxiliaires surf a- 
ciques dependant de parametres choisis dudit objet, stock^s. 

11. Dispositif selon I'une des revendications 1 Si 10, caractS- 
rise en ce que ledit module de forces internes (16) est propre 
& estimer les d^placements des noeuds du maillage volumique h 
partir du d6placement tir6 de 1' action d6finie et de forces 
externes, en particulier de type force gravitationnelle. 

12. Dispositif selon I'une des revendications 1 k 11, caract6- 
rise en ce que lesdits dSplacements estimSs des noeuds, autres 
que ceux de ladite maille surfaci^e comprenant ledit point 
d' intersection, sont obtenus par une m^thode choisie parmi au 
moins la mSthode d'Euler et la mSthode de Runge*Kutta. 

13. Dispositif selon la revendication 12, caract6ris6 en ce que 
ladite deformation est obtenue par la m6thode de Runge-Kutta 
dite " d ' ordre 4 " . 

14. Dispositif selon I'une des revendications 1 St 13, caracte- 
ris6 en ce que ledit module de forces internes (16) est propre, 
apr^s determination des d6placements estimes des noeuds, h 
supprimer au moins une liaison entre des noeuds voisins ou une 
maille volumique en fonction d'un premier critfere, puis h mettre 
h. jour les tenseurs de noeud et les tenseurs de liaisons en 
fonction de la ou des liaisons supprimees, et enfin h recalculer 
les forces internes des noeuds d'au moins la premiere partie du 
maillage volumxque. 

15. Dispositif selon la revendication 14, caracterise en ce que 
le premier crit&re porte sur au moins un paraun&tre choisi parmi 
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au moins une information transmise par ladite interf ace-utilisa- 
teur ( 4 ) et relative au type d ' outil (6,7) manoeuvre , une 
variation de volume de la maille volumique comprenant ladite 
liaison ou ladite maille volumique h. supprimer, et une variation 
de longueur d'une liaison de la maille volumique comprenant 
ladite liaison ^ supprimer. 

16. Dispositif selon la revendication 15, caract6ris6 en ce que 
ledit module de forces internes (16) est propre, aprSs determi- 
nation des d^placements estim^s des noeuds, ^ supprimer un noeud 
en cas de detection de suppression de toutes les liaisons qui 
relient ledit noeud aux noeuds voisins ou en fonction dudit 
premier crit^re, puis a mettre h jour les tenseurs de noeud et 
les tenseurs de liaisons en fonction du noeud et des liaisons 
supprim^s , et enfin ^ recalculer les forces internes des noeuds 
d'au moins la premiere partie du maillage volumique. 

17. Dispositif selon la revendication 16, caract6ris6 en ce que 
ledit module de forces internes < 16) est propr^, cas <ie 
suppression de liaison et/ ou de noeud et avant mise a jour des 
tenseurs de liaison et de noeud, k a j outer de nouveaux noeuds et 
de nouvelles liaisons de manigire Sl reroailler ledit maillage 
volumique . 

18. Dispositif selon I'une des revendications 1 h 17, caract6- 
ris6 en ce que le module de forces internes (16) est agenc6 pour 
estimer la force interne exerc6e sur chaque noeud dudit maillage 
volumic[ue dudit objet. 

19. Dispositif selon I'une des revendications 1 St 17, caract6- 
ris4 en ce que le module de forces internes (16) est propre ^ 
determiner les forces internes exerc^es sur les noeuds d'au 
moins une seconde partie dudit maillage volumique h partir de 
conditions aux limites dSfinies par des noeuds dits de connexion 
places & 1' interface entre les premiere et seconde parties, et 
d'une table de tenseurs de deformation, dont chaque tenseur est 
representatif de 1' influence d'un deplacement ei6mentaire de 
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chaque noeud d'au moins ladite seconde partie sur chaq[ue autre 
noeud d'au moins cette seconde partie. 

20. Dispositif selon la revendication 19, caract6ris6 en ce que 
5 les conditions aux limites servant au calcul des forces internes 
de la seconde partie sont d^finies par les forces internes 
calcul^es pour lesdits noeuds de connexion en tant que noeuds de 
la premiere partie. 

10 21. Dispositif selon la revendication 20, caractSris^ en ce que 
le module de forces internes (16) est agenc4 pour dSduire des 
valeurs des forces internes exerc^es sur les noeuds de la 
seconde partie du maillage volumique des valeurs de dSplacement 
des noeuds de connexion de maniere ci fournir des conditions aux 

15 limites permettant ci leur tour de calculer les forces internes 
des noeuds de la premifere partie. 

22. Dispositif selon la revendication 21, caractSris^ en ce que 
le module de forces internes (16) est propxe &: recalculex les 

20 valeurs des forces internes exerc^es sur les noeuds des premiere 
et seconde parties du maillage voluioique h partir des conditions 
aux limites d6termin6es success ivement sur les noeuds de 
connexion, jusqu'i obtention d'une position dite "d'6quilibre" 
des forces internes desdits noeuds de connexion. 

25 

23. Dispositif selon I'une des revendications 19 ^ 22, caract6- 
ris6 en ce que les parties du maillage volumique sont d^termi- 
n6es ci partir d'un critfere pr6d6termin6 portant au moins sur un 
param^tre des donndes-image dudit objet choisi parmi des 

30 param^tres physiques et des param^tres anatomiques, en particu- 
lier sur une intensity. 

24. Dispositif selon I'une des revendications 19 Si 23, caract6- 
ris6 en ce que la seconde partie est la partie compl6mentaire de 

35 ladite premiere partie. 

25. Dispositif selon I'une des revendications 19 d 24, caract^- 
ris6 en ce qu'il comprend un module de partitionnement () propre 
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h subdiviser ledit maillage volumique en lesdites parties, dont 
lesdites premiere et seconde parties. 

26. Dispositif selon I'une des revendications 1 ci 25, caract6- 
5 ris6 en ce que ledit module de determination de "collision" (18) 

est propre ^ determiner une collision entre au moins deux outils 
g6r6s par ladite interf ace-utilisateur (4). 

27. Dispositif selon la revendication 26, caract6ris6 en ce que 
10 chaque outil comprend au moins une extr6mit6 d' interaction (7) 

avec I'objet, laquelle est representee par au moins un point. 

28. Dispositif selon la revendication 27, caract6ris6 en ce que 
pour la tige de chaque outil est representee par une multipli- 

15 cite de points reliSs entre eux, ainsi qu'Si la ou les extr^mites 
( 7 ) , par des segments . 

29. Dispositif selon I'une des revendications 27 et 28, 
caracterise en oe que ledit module de determiTtation de collisi-on 

20 (18) est propre a cr^er un espace tridimensionnel englobant 
I'enveloppe externe dudit objet, puis h decomposer ledit espace 
en cellules volumiques dont le nombre est choisi de sorte que 
chaque cellule comporte un nombre de noeud dudit maillage 
volumique de 1' objet sensiblement egal au nombre de noeuds des 

25 autres cellules, chaque cellule intersectant ladite surface 
externe comport ant au moins un noeud, et en fin Si stocker dans 
des multiplets les coordonnees de chaque noeud en reference St la 
cellule volumique qui I'englobe. 

30 30. Dispositif selon la revendication 29, caracterise en ce que 
ledit module de detection de collision (18) est propre h 
determiner la presence d'un point de 1' outil dans ledit espace 
par une comparaison entre lesdits multiplets et les coordonnees 
dudit point, puis la cellule volumique dudit espace dans 

35 laquelle se trouve ledit point. 

31. Dispositif selon la revendication 30, caract6rise en ce que 
ledit module de detection de collision (18) est propre & 
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determiner la distance qui s6pare ledit point dudit ou desdits 
noeuds englob6s dans ladite cellule volumique de fagon a 
determiner la plus petite de ces distances, appeiee distance 
minimale, puis a determiner la distance qui s^pare ledit point 
5 dudit ou desdits noeuds englobSs dans un nombre predetermine de 
cellules volumiques voisines de la cellule volxamique dans 
laquelle il se trouve de f agon ci comparer ses distances h ladite 
distance minimale, puis determiner 1' ensemble des mailles 
surfaciques adjacentes au noeud associe h, ladite distance 

10 minimale de fagon h determiner si un segment defini par ladite 
position du point de I'outil et par sa position pr6c6dente 
intersecte I'une au moins de ces mailles surfaciques adjacentes, 
ce qui permet d'estimer les coordonnees (barycentriques ) du 
point de collision entre ledit objet et ledit outil, en vue de 

15 leur transmission audit module de forces internes (16). 

32. Dispositif selon la revendication 31, caracteris6 en ce que 
le module de detection de collision (18) est propre h modifier 
le contenu desdits multiplets entre deux determitiations de 

20 presence de points & l'int6rieur desdites cellules volumiques, 
en cas de detection d'une deformation par ledit module de forces 
internes ( 16) • 

33. Dispositif selon I'une des revendications 1 St 32, caracte- 
25 rise en ce que 1 ' interf ace-utilisateur (4) comprend un ensemble- 

harnachement (5) manoeuvrable par au moins une main d'op6rateur 
pour simuler la manoeuvre dudit outil. 

34. Dispositif selon I'une des revendications 1 Si 33, caracte- 
30 rise en ce qu'il comprend des moyens d'affichage (3) couples au 

module de raf raichissement (17) pour afficher en temps r6el sous 
forme d' image au moins les donn6es-image de 1' objet et d'une 
representation de 1' outil.. 

35 35. Dispositif selon I'une des revendications 1 & 34, caract6- 
rise en ce que le jeu de donnees-image represente une image 
numerique tridimensionnelle d'une region comportant au moins un 
objet 3D, dont 1' objet designe. 
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36- Dispositif selon la revendication 35 ^ caract6ris6 en ce que 
1' image num^rique est une image m6dicale. 

37- Proced^ ^lectronique de traitement de donn6es-image, du 
5 type comprenant les 6 tapes : 

- pr6voir une interf ace-utilisateur (4) capable de generer un 
retour de force, conf orm6ment aux reactions d'un outil (6,7), 

- ^tablir, h partir d'une loi de deformation et d'une action 
d6finie par 1 ' interf ace-utilisateur (4) et representative d'une 

10 manoeuvre dudit outil, un champ de forces internes entre des 
noeuds d'un maillage volumique fonction d'un maillage surfacique 
d'un objet 3D apparaissant dans un jeu de donnSes- image, 

- determiner la force de reaction de 1' objet qui correspond h sa 
deformation estimee Si partir des forces internes, de sorte que 

15 la force gen6ree par 1' interf ace-utilisateur (4) soit sensible- 
ment equilibree par cette force de reaction, 

- calculer de nouvelles donnees-image de 1' objet, en presence 
des deformations estimees , compietees de la representation dudit 
outil , 

20 caracterise en ce que I'on prevoit une 6tape dans laquelle 

on estime un point d' intersection entre une droite mater ialisant 
un deplacement tire de ladite action def inie et ledit maillage 
surfacique, et en ce que I'on estime la force interne exercee 
sur les noeuds d'une premiere partie au moins du maillage 

25 volimique de 1' objet, h partir du deplacement applique aux 
noeuds appartenant a la maille surfacique contenant ledit point 
d' intersection, de conditions aux limites, et de tenseurs de 
noeud et de tenseurs de liaison issus respectivement pour chaque 
noeud et chaque liaison de cette partie au moins, de matrices de 

30 rigidite propres h chaque maille volumique d'au moins ladite 
premiere partie et fonction de la loi de deformation • 



BNSDCXID: <FR 2771202A1J_> 



2771202 



1/7 




//////// 



.r-r. 



FIG.1 



LIE 




MEMOIRE 
PARTAGEE 



I 



RESEAU 



POSITIONS 
(x.y.z) 



CALCUL 
FORCES 
INTERNES 



T 



I 
I 

I 
I 



RESEAU 



I 



CALCUL 
DEFORMATION 



FORCES 
(F^,Fy.Fz) 



POSITIONS 
(x.y.z) 



18 

rrr- — 

jj DETECTION DE 
COLUSIONS 



BUS ISA 



I 



FIG.2 




STATION GRAPHIQUE 



UNITE DECONTROLE 



BNSCXXJID: <FR 2771202A1_L> 




BNSOOCID: <FR 2771202A1_I_> 





F1G.9 



11122233ai4i55S666777ft8S999AAABBBCCCDDDC.tLFFr 

1XXXX..X.X. ..xxxxx xxxxxx 

iXxx.xx.x. . . .xxxxx xxxxxx 

ixxx.xxx«x...xxx..x xxxxxx 

2X. .XXX. . .X.XXXX.. .XX XXX .. .XXX 

2 « JC X ^ IC m m • 9 ^ ^ • • • ^ 5^ ^ « * « ^ "JCX^^ • • 

. .X.XXXX X XXX. ..XXX 

Sx.x.. .XXXX..XXX XX XXX.. .xxxxxx 

3.x. . ..xxx.xxxxx XX XXX.. .xxxxxx 

3X.X. . .xxx.xxxxx X XXX. . .xxxxxx 

a.X.XX. .xxxxxx. • • vXX. ..XXX«. • .XXX •« 

i.. - .X. .XXXXXXXX XX. . , .XXX XXX , 

4...X.X.XXXXXXXX X .. .XXX XXX. . . 

5XXXXXXXXXXXXXXX xxxxxxxxxxxxxxx 

sxxxxxxxxxxxxxxx xxxxxxxxxxxxxxx 

5XXXXXXXXXXXXXXX xxxxxxxxxxxxxxx 

eXX XXXX. .X.X. ..XXX. . .xxxxxx.. .XXX 

6XX XXX.XX.X. .XXX. ..xxxxxx. . .XXX 

S. .X xxx.xxx.x.. .XXX. ..xxxxxx.. .XXX 

7.. .XX X..XXX...X.XXXXXXX...XXX 

7.. .XX XXXXX....X. xxxxxx. ..XXX 

7 X XXXXX...X.XXXXXXX.. .XXX 

I XX X.X. ..XXXX.. XXX.. .XXX.. .xxxxxx 

M f . . . XX. a . . . . .X.. • .xxx.xxxxx^. .XXX .xxxxxx 

« X X.X. ..xxx.xxxxx. XXX.. .xxxxxx 

9 XX X.XX..XXXXXXXXX xxxxxx 

0 XX X..XXXXXXXXXXX xxxxxx 

1 X X.X.XXXXX.XXKXM9II XXXXXX 

A xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

A xxxxxxxxxxxxxxxxxxxxxxxxxxxxxs 

A xxxxxxxxxxxxxxx&xxxxxxxxxxxxxx 

e. . . XXX. . .xxxxxx. . . XXX. . .xxxxxxxxx. . .xxxxxxxxs 

B. ..XXX... xxxxxx... XXX.. .XXXXXXXXX...XXXXXXXXX 

?. ..XVX;. .XXX&XX...XXX.. .XXXXXXXXX...XXXXXXXXX 
XXX. . .XXX. . .xxxxxx. . .XXX. . .XXX.. .XXKXXXXXXXXX 
CXXX. . . XXX. . .xxxxxx. . . XXX. . .XXX. . .xxxxxxxxxxxx 

SXXX. . .XXX. . . xxxxxx. . .XXX. . . XXX. . .XX.KXXXXXXXXX 
xxxxxx XXXXXXXXX xxxxxxxxxxxx. . . XXX 

B xxxxxx XXXXXXXXX XXKXXXXXXXXX. ..XXX 
xxxxxx XXXXXXXXX xxxxxxxxxxxx. ..XXX 

t XXXXXXXXX xxxxxxjrtrxxvJrjTArjv'. . .xxxxxx 

r XXXXXXXXX xxx'xxxxxxxx'xxxx. ..xxxxxx 

r XXXXXXXXX xxxxxxxxxxxxxary. . .xxxxxx 

r XXX. . .xxxxxx. ..xxxxxxxxxxxxxxxxxxxxxxxx 

r XXX.. .xxxxxx... xxxxxxxxxsrxxxxxxxxxxxxxx 

r XXX. . .xxxxxx. . .xxxxxxxxxxxxxxxxxxxxxxxx 



BNSDOCID: <FR 277l2QgA1 i > 



> 

2771202 





4/7 


MAILLAGE 




SURFACIQUE 








MAILLAGE 


OEFlNfTlON DESUTRAEDR 


VOLUMIQUE 


ET DES TRIANGLES 



1 

MATRICES DE 
RIGIDITE 
ELEMENTAIRES 



I 

MATRICE DE 
RIGIDITE 

GLOBALEIK] 



MAILLAGE FINAL 
(MASSESJRIANGLES, 

TETRAEDRES)+(LIAlSONS) 



DEFINITION 
DES LIAISONS 



FIG.10 



CALCUL DEL ACCELERATION: 
SOMME DES FORCES 



( integration) 



miseaxxjr 
de la scene 



CALCUL DE LA VITESSE 



(^INTEGRATIONJ 



CALCUL DE LA 
NOUVELLE POSITION 



miseajouroestenseursJ 









1 oufy 



NON 

NOMBRE 
D'ITERATOMS 

>SEUIL 



JEMPSDECALCUL' 
>DELAIMAXI 

V 




SUPPRESSION 
D'UNE MASSE 



SUPPRESSION 
D'UNE LIAISON 



RG.11 



BNSDOCID: <FR_2771202A1 _!_> 




CDL= FORCES AUX 
NCEUDSDE CONNEXION 



INTERAQIONS 



I 




CDL=DEPLACEMENTSDtS 
NCEUDSDE CONNEXION 



MAILLAGE Al 



I 



INTERACTIONS 



FIG.15 




BNSDCX:iD: <FR . 2771202A1J_> 



2771202 



REPUBUQUE FRANQAISE 



INSTTTUT NATIONAL 
dela 

PROPRIETE INDUSTRIELLE 



RAPPORT DE RECHERCHE 

PRELIMINAIRE 

etabll 8ur la base dee demieres revendicatione 
depoeees avant le commenoement de la recherche 



rr d'enregistremant 
national 



FA 552415 
FR 9714506 



Categoric 



DOCUMENTS CONSIDERES COMME PERTINENTS 



Citation du dooument avao indication, an c 
des parlias perttnentaa 



( de beaoin. 



Revendi&ations 
oonceTn6ee 
de la demande 
examine 



WO 97 19440 A (IMMERSION HUMAN INTERFACE 
CORP) 29 mai 1997 

* le document en entier * 

WO 96 30885 A (GILLIO ROBERT G) 3 octobre 
1996 

* le document en entier * 

GB 2 298 931 A (MARCONI GEC LTD) 18 
septembre 1996 

* le document en entier * 

US 5 625 576 A (MASSIE THOMAS H ET AL) 29 

avril 1997 

* le document en entier * 

WO 96 16389 A (STANEFF JOHN E JR ;MOORE 
ROBERT S (US); HARTHAN LEWIS JOHN III (US) 
30 mai 1996 

* le document en entier * 



1-37 



1-37 



1.37 



1.37 



1.37 



DOMAINES TECHNIQUES 
RECHERCHES (lnt.CL.6) 



G05B 

B25J 
G06F 
G09B 



Dale d'acMvement de la rechecohe 

10 aout 1998 



Hauser, L 



CATEGORIE DES DCXiUMEhTTS CfTES 

X : parliouK^mant pertinant A tui aaul 

Y : par1iouli6raman1 pwtinant an oembfnaiaon avM un 

autre dooumant da la nn*ma catAgoria 
A : partinant k Tanoontra d'au moina una ravandioation 

ou axn^f^jHan taohnolo0iqua gAnaral 
O : divulgation non-Aonta 

P : doownant intareataira 



T : th^rie ou prinoipa h la baae da Tinvanbon 

E : dooumant da brevet banafioiant d'una data ant6neura 

. A ta data de d6pOt at qui n'a 6ta pubM qu'& oette data 

de d6pdt ou qu'A une data poat^haura. 
0 : oit^ dana ta demand* 
L : oct6 pour (f autres raiaona 

a : mafrbia da la rnftma famila, dooumant oorraipondant 



SNSDOCID: <FR 2771202A1_L> 



